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This article presents a review of the field of inclusive quasi-elastic electron-nucleus scattering. It 
discusses the approach used to measure the data and includes a compilation of data available in 
numerical form. The theoretical approaches used to interpret the data are presented. A number 
of results obtained from the comparison between experiment and calculation are then reviewed. 
The analogies and differences to other fields of physics exploiting quasi-elastic scattering from 
composite systems are pointed out. 
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I. INTRODUCTION 

The energy spectrum of high-energy leptons (electrons 
in particular) scattered from a nuclear target displays a 
number of features. At low energy loss (u>) peaks due 
to elastic scattering and inelastic excitation of discrete 
nuclear states appear; a measurement of the correspond- 
ing form factors as a function of momentum transfer |q| 
gives access to the Fourier transform of nuclear (transi- 
tion) densities. At larger energy loss, a broad peak due 
to quasi-elastic electron-nucleon scattering appears; this 
peak — very wide due to nuclear Fermi motion — cor- 
responds to processes where the electron scatters from 
an individual, moving nucleoli, which, after interaction 
with other nucleons, is ejected from the target. At even 
larger to peaks that correspond to excitation of the nu- 
clcon to distinct resonances are visible. At very large 
u), a structureless continuum due to Deep Inelastic Scat- 
tering (DIS) on quarks bound in nucleons appears. A 
schematic spectrum is shown in Fig. [1] At momentum 
transfers above approximately 500 MeV/c the dominant 
feature of the spectrum is the quasi-elastic peak. 

A number of questions have been investigated using 
quasi-elastic scattering: 

• The quasi-elastic cross section integrated over 
electron energy loss is proportional to the sum 
of electron-nucleon cross sections. Historically, 
this has been exploited in order to measure the 
neutron charge and magnetic form factors using 
mainly light (A < 4) nuclear targets. Today the 
emphasis has shifted to exposing possible medium 
modifications of the nucleoli form factors. 
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FIG. 1 Schematic representation of inclusive cross section as 
function of energy loss. 



Another integral property of the quasi-elastic peak, 
its width, provides a direct measure of the average 
momentum of nucleons in nuclei, and has been used 
to determine nuclear Fermi momenta; contrary to 
other observables such as densities, quasi-clastic 
scattering provides a direct determination via 
an observable sensitive to the momenta of nucleons. 



The shape of the quasi-elastic peak depends on the 
distribution in energy E and momentum k of the 
initially bound nucleons. Precise measurements 
give indirect access to (integrated properties of) 
the nuclear spectral function S{\a, E) describing 
this distribution. In particular, the tail of the 
quasi-elastic peak at large |q| and low u is 
sensitive to the tail of the spectral function at 
large |k|. 



• Processes more complicated than one-nucleon 
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knockout, in particular those related to non- 
nucleonic degrees of freedom and Meson Exchange 
Currents (MEC), also play a role. They can be in- 
vestigated by separating the quasi-elastic response 
into the longitudinal (charge) and transverse 
(magnetic plus convection) pieces, the latter being 
preferentially affected by MEC. 

• Scaling is one of the unique features of quasi-elastic 
scattering. This refers to the fact that the inclusive 
cross section, which a priori is a function of two in- 
dependent variables, |q| and uj. depends on a sin- 
gle variable y(q,u>). This scaling property, a conse- 
quence of the kinematics of the underlying electron- 
nucleon elastic scattering process, provides a strong 
handle on the reaction mechanism. Further, the 
scaling violations that are observed reveal how the 
dynamics go beyond the Impulse Approximation 
(IA) picture of quasi-elastic scattering. 

Inclusive quasi-elastic electron-nucleus scattering is 
not the only process of this type. This quasi-elastic pro- 
cess occurs in various other areas of physics, and is being 
exploited to learn about the dynamics of the underlying 
composite system: 

- Quasi-elastic scattering of keV-energy photons from 
electrons bound in atoms provides information about the 
energy and momentum distribution of bound electrons. 

- Scattering of eV to kcV neutrons from condensed- 
matter systems such as liquid Helium provides a mea- 
surement of the He momentum distribution and correla- 
tions, and has been exploited for a long time to isolate 
the effects of the Bose condensate in superfluid Helium. 

- Quasi-elastic scattering of GeV-energy leptons on the 
quarks bound in nucleons (DIS) has provided a wealth of 
information on the quark distribution functions; the ob- 
servation of scaling violations in DIS has taught us much 
about the dynamics of strong interactions. 

Inclusive electron scattering from nuclei is a subclass of 
quasi-elastic processes, the most obvious other represen- 
tative being (e, e'p). When compared to exclusive pro- 
cesses where the knocked-out nucleon is detected, (e, e') 
corresponds to an integral over all final states of the nu- 
cleon, and consequently provides less specific informa- 
tion. On the other hand, (e, e') is more directly related 
to the dynamics of the initial hadronic (nuclear) ground 
state. The complications of the final hadronic contin- 
uum play a much smaller role. This is true particularly 
at large |q|, as the electron is sensitive to the fate of 
the recoiling nucleon in a region of size l/|q|; the Final 
State Interaction (FSI) beyond that region affects only 
the more exclusive processes. 

Quasi-elastic scattering from nuclei, compared to the 
other quasi-elastic processes mentioned above, has one 
distinct drawback one must deal with: the nucleon the 
electron scatters from is not elementary, but can be ex- 
cited to various states. At large momentum transfer the 
result is that only the low-w side of the quasi-elastic 



peak can be exploited, the large-o; side is obscured by 
the overlap with A excitation. 

During the last decade, the remarkable progress of our 
understanding of quasi-elastic scattering has allowed us 
to define a number of features precisely. In particular 
for light nuclei and nuclear matter, quasi-elastic scatter- 
ing has provided accurate determinations of (integrated 
properties of) the nuclear spectral function. Detailed in- 
vestigations of the reaction mechanism and the condi- 
tions necessary/sufficient for scaling have clearly shown 
that the dominant process is elastic scattering from indi- 
vidual nucleons. The different scale-breaking mechanism 
have been identified. The effects beyond IA, mainly due 
to FSI, have come under much better control. The kine- 
matical region where more complicated processes, such as 
MEC, are important has also been identified and the size 
of two-body current contributions has been understood. 

In this review we try to give a reasonably compre- 
hensive discussion of the various aspects of quasi-elastic 
scattering. The quantitative understanding of the cross 
section starts from a description in terms of IA; effects 
beyond IA such as the role of FSI and MEC are addressed 
next. We then describe some of the experimental aspects 
important for study of this reaction. We then give a com- 
pilation of the experiments that have been performed and 
point out where the cross sections (most often not pub- 
lished in numerical form) can be found. We also briefly 
describe experiments that provide the nucleon form fac- 
tors needed to understand quasi-elastic scattering. We 
subsequently discuss scaling and the related superscal- 
ing. For light nuclei and nonrelativistic final states, exact 
calculations can be performed. For the lower momen- 
tum transfers, an alternative approach, the use of the 
Euclidean response, is available and presented. We then 
study the results obtained after a longitudinal/transverse 
(L/T) separation of the cross section, and their impact 
on the Coulomb sum rule. A bothersome correction, the 
effect of Coulomb distortion on the cross sections, is ad- 
dressed as well. We also show how data for an important 
model system for nuclear theory, infinite nuclear matter, 
can be obtained. Last, we address other fields of quasi- 
elastic scattering, and discuss the common aspects. 



II. ELECTRON-NUCLEUS SCATTERING IN IMPULSE 
APPROXIMATION 

A. Electron-nucleus cross section 

The differential cross section of the process 

e + A^e'+X , (1) 

in which an electron of initial four-momentum k e = 
(E e ,\s. e ) scatters off a nuclear target to a state of four- 
momentum k' e = (E e i , k e ' ), the target final state being 
undetected, can be written in Born approximation as 
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(|ltzvkson and Zubeij . [l980h 

d 2 a _ a 2 E e > 
dtt e ,dE e , ~ Q 1 



(2) 



where a = 1/137 is the fine structure constant, dfl e i is 
the differential solid angle in the direction specified by 
k e ', Q 2 = —q 2 and q = k e — k e i = (w,q) is the four 
momentum transfer. 

The tensor L^, that can be written neglecting the 
lepton mass as 



L 



1 1 V 



2[k^k u e ,+k v e k^-g^(k e k e ')] , 



(3) 



where g^ v = diag(l, — 1, — 1, — 1) and (k e k e >) = E e E e i — 
k e ■ k e < is fully specified by the measured electron kinc- 
matical variables. All the information on target structure 
is contained in the tensor W^, whose definition involves 
the initial and final nuclear states |0) and \X), carrying 
four-momenta po and px , as well as the nuclear current 
operator J M : 



E 

X 



(0\J»\X}{X\r\0)SW(p o + q- 



Px, 



(4) 



where the sum includes all hadronic final states. 

The most general expression of the target tensor of 
Eq. (dJ), fulfilling the requirements of Lorentz covariance, 
conservation of parity and gauge invariance, can be writ- 
ten in terms of two structure functions W\ and Wi as 



Wi 
W2 



Oog ) 

q 2 



-q' 



Po 



(poq) 



q 



-q 



(5) 



where m is the nucleon mass and the structure functions 
depend on the two scalars Q 2 and (poq). In the target 
rest frame (poq) — muj and W\ and Wi become functions 
of the measured momentum and energy transfer |q| and 

UJ. 

Substitution of Eq. ([5]) into Eq. ^ leads to 
da 



dQ„i dE P 



dflp 



M 



M/ 2 (|q|, W ) + 2VF 1 (|q|,o;)tan 2 - 



(6) 



where 6 and (da/dil e ,) M = a 2 cos 2 (6>/2)/4£ e sin 4 (6»/2) 
denote the electron scattering angle and the Mott cross 
section, respectively. 

The right hand side of Eq. ^ can be further rewrit- 
ten singling out the contributions of scattering processes 
induced by longitudinally (L) and transversely (T) po- 
larized virtual photons. The resulting expression is 



dQ,pi dEp 



da 
dflpi 



M 



Q4 

rjj RU\q\,v) 



(7) 



where the longitudinal and transverse structure functions 
are trivially related to W\ and Wi through 



i?r(|q|,w) = 2Wi(|q| 



(8) 



and 



^lfl L (|q|,a,) - W 2 (\q\,u) - j?-Wi(\q\,u) . (9) 

In principle, calculations of W^ v of Eq. ([4| at mod- 
erate momentum transfer (|q| < 0.5GeV/c) can be car- 
ried out within nuclear many-body theory (NMBT), us- 
ing nonrelativistic wave functions to describe the initial 
and final states and expanding the current o perator in 
powers of |q|/m (jCarlson and Schiavillal 1 19981) . where m 
is the nucleon mass. The available results for medium- 
heavy targets have been mostly obtained using the mean 
field approach, supplemented by the inclusion of model 
residual i nteractions to take into account long range cor- 
relations ( Dellafiore et q/J . Il985f ). 

On the other hand, at higher values of |q|, correspond- 
ing to beam energies larger than ~1 GeV, describing the 
final states \X) in terms of nonrelativistic nuclcons is no 
longer possible. Due to the prohibitive difficulties in- 
volved in a fully consistent treatment of the relativistic 
nuclear many-body problem, calculations of W^ v in this 
regime require a set of simplifying assumptions, allowing 
one to take into account the relativistic motion of final 
state particles carrying momenta ~ q, as well as inelas- 
tic processes leading to the production of hadrons other 
than protons and neutrons. 

The IA scheme is based on the assumptions that i) as 
the space resolution of the electron probe is ~ l/|q|, at 
large momentum transfer scattering off a nuclear target 
reduces to the incoherent sum of elementary scattering 
processes involving individual bound nucleons (see Fig. [5]) 
and ii) there is no FSI between the hadrons produced in 
electron- nucleon scattering and the recoiling nucleus . 
Under these assumptions, a relativistic particle in the fi- 
nal state is completely decoupled from the spectator sys- 
tem, and the description of its motion becomes a simple 
kincmatical problem. 

Within the IA picture the nuclear current appearing 
in Eq. ^ is written as a sum of one-body currents 



(10) 



while I A) reduces to the direct product of the hadronic 
state produced at the electromagnetic vertex, carrying 
momentum p^, and the state describing the (A — 1)- 
nucleon residual system, carrying momentum = q — 



1 Coherent effects, not included in the impulse approximation 
picture, appear in DIS even at large |q| for values of the Bjorkcn 
scaling variable x = Q 2 /2mui £ 0.2, corresponding to very large 
electron energy loss. 
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FIG. 2 Schematic representation of the IA regime, in which 
the nuclear cross section is replaced by the incoherent sum 
of cross sections describing scattering off individual nucleons, 
the recoiling (A — l)-nucleon system acting as a spectator. 



p x (in order to simplify the notation, spin indices will be 
omitted) 



\X) -> \x,p x }® \K,p n ) 
Using Eq. (fTT|) we can replace 



(11) 



^2\X)(X\ [ d 3 Px \x,p x )(p x ,x\ 

X x J 

X ^d 3 p K |^,PK)(P7Z,^| . (12) 



R 



Substitution of Eqs. fT0 j) - ([T2"]) into Eq. Q and insertion 
of a complete set of free nucleon states, satisfying 



J d 3 k |N,k)(k,N| = 1 , 



(13) 



then leads to the factorization of the nuclear current ma- 
trix clement according to 

/ \ V2 



(0\J»\X) 



\/\pn\ 2 + m 2 



<0|7e,p w ;N,-p w ) 



Y,(-VK,N\tf\x,p x ) 



(14) 



where the factor (raj \J\pn\ 2 + m 2 ) 1 ' 2 takes into account 
the implicit covariant normalization of (—pn,N\ in the 
matrix element of . 

As a result, the incoherent contribution to Eq. (j4|) can 
be rewritten in the form 

W» v = [ d 3 Pn d 3 Px \(0\ll,PK;N,-p n }\ 2 

x I](-p^ N bTI^ p*)(p*. z|#|n,-ptc) 

i 

x 5^{q-p K -p x )5{cj + E -E n -E x ), (15) 

where Eq is the target ground state energy and E-r. = 
\/\piz\ 2 + M n , M-ji being the mass of the recoiling sys- 
tem and E x the energy of the final state X. 
Finally, using the identity 

<S(w + E -E n - E x ) = J dE S(E -m + E - E n ) 
x 5(uj-E + m-E x ) , (16) 



and defining the nucleon spectral function as 
S N (k 7 E) = ^|(0|ft,-k;N,k}| 2 



■R 



x 8{E -m + E - E n ) , 



(17) 



where the index N — p,n labels either a proton or a 
neutron, we can cast Eq. (U) in the form 



W^(q,uj) = J d 3 kdE f^j [ZS p (k,E)w£ v 

+ (A-Z)S n (k,E)w£ u ] , (18) 

A and Z being the target mass number and number of 
protons, respectively. In Eq. (fl8|) . E^ = \Z\k 2 \ + m 2 and 

w% = ^(k,N|^|x,k + q)(k + q,.x| J X r |N,k) 

X 

x 5(u + E k -E x ) . (19) 
with (sec Eqs. (JTSj) and (fTTjl) 

2 = E x -E u = E Q +uj-E n -Ek = oj-E+m-E^ . (20) 

The above equations show that within the IA scheme the 
definition of the electron-nucleus cross section involves 
two important elements: i) the tensor , defined by 
Eq. (fl"9|) . describing the electromagnetic interactions of 
a bound nucleon carrying momentum k and ii) the spec- 
tral function, defined by Eq. (|17|) , yielding its momentum 
and removal energy distribution. These quantities will be 
further discussed in the following sections. 



B. Electron scattering off a bound nucleon 

While in electron-nucleon scattering in free space the 
struck particle is given the entire four momentum transfer 
q = (u>, q), in a scattering process involving a bound 
nucleon a fraction Slu of the energy loss goes into the 
spectator system. This mechanism emerges in a most 
natural fashion from the IA formalism. 

Assuming that the current operators are not modi- 
fied by the nuclear environment, the quantity defined 
by Eq. (fl"9)) can be identified with the tensor describing 
electron scattering off a free nucleon at four momentum 
transfer q = (q, ui). Hence, Eq. (fT9|) shows that within IA 
binding is taken into account through the replacement 



q = (w, q) -> q = (w, q) 



(21) 



The interpretation of dcu = lu — uj as the amount of 
energy going into the recoiling spectator system becomes 
particularly transparent in the limit |k|/m <C 1, in which 
Eq. $2Q§ yields Slo = E. 

In the case of quasielastic scattering, to be discussed 
in this review, of Eq. (|19[) can be obtained from the 
general expression (compare to Eq. 



w 



N 



w 



N 



/??- 



-9 



(kq). 

q 2 



k v - 



q 2 



(22) 
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where k = (i?k,k) and the two structure functions W\ 
and u>2 are simply related to the measured electric and 
magnetic nucleon form factors, Ge n and Gm n , through 



,.N 



,.N 



4m 2 \ 2m 



1 - q 2 /4m 2 



G 2 



Ma 



1_ 

2m 



(23) 



Gi., - 



4m 2 



G 2 



Ma 



(24) 



While the replacement of to with u) is reasonable on 
physics grounds and in fact quite natural in the context 
of the PWIA analysis presented above, it poses a consid- 
erable conceptual problem in that it leads to a violation 
of current conservation, which requires 



. 



(25) 



A widely used prescription to overco me this di f ficulty has 
been proposed in the early 1980's ( de Forestl Il983f ). It 



amounts to using a tensor whose time components 



are given by Eq. (|19p . while the longitudinal ones are 
obtained from the time components requiring that the 
continuity equation (|25[) be satisfied. Taking the z-axis 
along the direction of q one obtains 



w 



(q) , for fi and/or v = 



N 



(<7) 



(26) 



with q given by Eq. (|2"Tj) . 

The above prescription is manifestly non-unique. For 
example, one might have just as well chosen to use 
Eq. (fl"9)) to obtain the longitudinal components of w^ v 
and the continuity equation to obtain the time compo- 
nents. However, differences between these two proce- 
dures to restore gauge invariance only affect the longi- 
tudinal response Rl- As a consequence, they are ex- 
pected to become less and less important as the momen- 
tum transfer increases, electron scattering at large |q| be- 
ing largely dominated by transverse contributions. Nu- 
merical calculations of the quasielastic cross section off 
3 H e have confirmed that these differences become indeed 
vanishingly small already i n the pion production region 
(|Meier-Haiduk et aZ.I . [l989h . 

The uncertainty associated with the implementation of 
current conservation has been recently analyzed by Kim, 
who has compared the inclusive cross sections computed 
using the Coulomb gauge (equivalent to the deForest pre- 
scription) to those obtained from the Landau and Weyl 
gauges for a variety of targets and kincmatical condition 
(|Kiml . l2004h . The results corresponding to the different 
choices agree with one another to better than ~ 1 % 
for backward angle kinematics, thus confirming that the 
transverse response is unaffected by gauge ambiguities. 
For the longitudinal response, Kim finds more substantial 



ambiguities, as a consequence of the huge scalar/vector- 
potentials present in the a — lu model; these large poten- 
tials, however, only occur when trying to describe nuclei 
— strongly correlated systems — using a mean-field ba- 
sis. 

Other prescriptions bas ed on the use of free nucleon 
cur rents (|de Forestl . Il983f) have been recently employed 
in the analysis of (e, e'p) data in the region of large miss- 
ing momentum and energy, where the effects of using dif- 
ferent off-shell extra polations of electron-nucleon scatter- 
ing become sizable (|Rohd . |2004[) . More fundamental ap- 
proaches, involving the explicit calculation of the off-shell 
form fac tors, necessarily rely on oversimplified dy namical 
models (|Naus and KochL Il987t iNaus et all Il990h . 

In conclusion the violation of gauge invariance in the 
IA scheme, while in principle an intricate one, turns out 
to be only marginally relevant to inclusive electron scat- 
tering at large momentum transfer. The main effect of 
nuclear binding can be easily accounted for with the re- 
placement LU — ► LU. 



C. The nuclear spectral function 

The spectral function, defined by Eq. (JT7J) , gives the 
probability of removing a nucleon from the target nucleus 
leaving the residual system with energy En = Eo—m+E. 

Within the shell model picture, based on the assump- 
tion that nucleons in a nucleus behave as independent 
particles moving in a mean field, the spectral function 
reduces to 



(27) 



e{F} 



where <p n (k) is the momentum-space wave function asso- 
ciated with the single particle shell model state n, E n is 
the corresponding energy eigenvalue and the sum is ex- 
tended to all occupied states belonging to the Fermi sea 
{F}- 

The results of electron- and hadron-induced nucleon 
knock-out experiments have provided overwhelming ev- 
idence of the inadequacy of the independent particle 
model to describe the full complexity of nuclear dynam- 
ics. While the peaks corresponding to knock-out from 
shell model orbits can be clearly identified in the mea- 
sured energy spectra, the corresponding strengths turn 
out to be consistently and sizably lower than expected, 
independent of the nuclear mass number. 

This discrepancy is mainly due to the effect of dynam- 
ical correlations induced by the nucleon-nucleon (NN) 
force, whose effect is not taken into account in the in- 
dependent particle model. Correlations give rise to scat- 
tering processes, leading to the virtual excitation of the 
participating nucleons to states of energy larger than the 
Fermi energy, thus depleting the shell model states within 
the Fermi sea. As a result, the spectral function acquires 
tails extending to the region of large energy and momen- 
tum, where SsAf(k, E) of Eq. (f2"T)) vanishes. 
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The typical energy scale associated with NN correla- 
tions can be estimated considering a pair of correlated 
nucleons carrying momenta ki and k2 much larger than 
the Fermi momentum (~ 250 MeV/c). In the nucleus 
rest frame, where the remaining A — 2 particles carry low 
momenta, ki « — k2 = k. Hence, knock-out of a nucleon 
of large momentum k leaves the residual system with a 
particle in the continuum and requires an energy 



thr 



2m 



(28) 



much larger than the typical energies of shell model states 
(~ 30 MeV). The above equation, where E t } lr denotes the 
threshold for two-nucleon removal, shows that large re- 
moval energy and large nucleon momentum are strongly 
correlated. 

Realistic theoretical calculations of the spectral func- 
tion have been carried out within NMBT, according to 
which the nucleus consists of a collection of A nucle- 
ons whose dynamics are described by the nonrelativistic 
Hamiltonian 



i=l 



2£ 

2m 



A 
j>i=l 



k>j>i=l 



(29) 



In the above equation, ki is the momentum of the i-th 
constituent and and Vijk describe two- and three- 
nucleon interactions, respectively. The two-nucleon po- 
tential, that reduces to the Yukawa one-pion-exchangc 
potential at large internucleon distance, is obtained 
from an accurate fit to the available data on the two- 
nucleon system, i.e. deuteron pro p erties and ~ 4000 
NN scattering data (jWiringa et all Il995h . The addi- 
tional three-body term Vijk has to be included in or- 
der to account for the binding energ ie s of the three- 



nucleon bound states ( Pudliner et M Il995h and the 



emp irical saturation properties of uni form nuclear mat- 
ter ( Akmal and Pandharipande! Il99?t ); this term results 
from the fact that non-nucleonic constituents (such as 
A's) have been excluded. 

The many-body Schrodinger equation associated with 
the Hamiltonian of Eq. (|2"9")l can be solved exactly, 
using stochastic methods, for nuclei with mass num- 
ber A < 12. The resulting energies of the ground 
and low-lying excited states are in excellent agree- 
ment with the experimental data ( Pieper and Wiringal . 
l200ll ). Accurate calculations can also be carried out 
for uniform nuclear matter, exploiting translational in- 
variance and using either a variational approach based 
on cluster expansion and chain s ummation techniques 
( Akmal and Pandharipande , 19971) . or G- matrix pertur- 



bation theory (jBaldo et al 



20001) . 



Nonrelativistic NMBT has been employed to ob 
tain the spectral functions of the three-nucleon sys 



terns (ICiofi degli Atti et adll980tlDieperink et al 



Meier- Ha iduk et al. , 19831 ). Oxygen ( Geurts et al 



1976; 



1996: 



Polls et all 119971) and symmetric nuclear matte r, hav 
ing A -> 00 and Z = A/2 (|Benhar et all Il989t 



iRamos et all Il989h . Calculations based on NMBT but 
involving some simpl ifying assumptions have been also 
carried out for 4 He ( Benhar and Pandharipande . 1993t 



ICiofi degli Atti et al.l . l!990t iMorita and Suzukil . Il99ll) . 



A S(k,E) 




180MeV 



FIG. 3 Nuclear matter spectral function ca lculated using 
corre lated basis function perturbation theory (jBenhar et all 
I1989I) . 

As an example, Fig. [3] shows the results of a nu- 
clear matter calculation 2 carried out using correlated ba- 
sis fu nction (CBF) perturbation theory ( Benhar et all 
1989). In addition to the peaks corresponding to sin- 
gle particle states, i.e. to bound one- hole states of the 
(A-l)-nucleon system, the resulting Sn(^,E) exhibits 
a broad background, extending up to E ~ 200 MeV 
and |k| ~ 800 MeV/c, associated with n-hole (n-1)- 
particle (A-l)-nucleon states in which at least one nu- 
cleon is excited to the continuum. The correlation ridge 
at E ~ k 2 /2m (see Eq. ([2"8|)) is clearly visible. Note that, 
in absence of interactions, the surface shown in Fig. [3] 
collapses to a collection of (^-function peaks distributed 
along the line \E\ = k 2 /2m, with |k| < k F « 250 MeV/c. 

The proton spectral functions of nuclei with A > 4 have 
been m odeled using t he Lo cal Density Approximation 
(LDA) (jBenhar et all Il994l ) , in which the experimental 



2 As in symmetric nuclear matter S p (k, E) = S re (k, E), the spec- 
tral function shown in the figure corresponds to an isoscalar nu- 
cleon. 
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information obtained from nuclcon knock-out measure- 
ments is combined with the results of theoretical calcula- 
tions of the nuclear matter S(k, E) at different densities. 

The kincmatical region corresponding to low miss- 
ing energy and momentum, where shell model dynamics 
dominates, has been extensively studied by coincidence 
(e, e'p) experiments. The spectral function extracted 
from the data is usually written in the factorized form 
(compare to Eq. l[27])) 



Smf(K E) 



ne{F} 



t (li)\ 2 F n (E-E n ) 



(30) 



where the spectroscopic factor Z n < 1 and the func- 
tion F n (E — E n ), describing the energy width of the n- 
th state, account for the effects of residual interactions 
not included in the mean field picture. In the Z n — ► 1 
and F n (E - E n ) -> 5(E - E n ) limit Eq. (30]) reduces to 
Eq. (27]). 

The correlation contribution to the nuclear matter 
spectral function has been calculated using CBF per- 
turbation the o ry fo r a wide range of density values 
(Benhar et all 11994] ). Within the LDA scheme, these 
results can be used to obtain the corresponding quantity 
for a finite nucleus of mass number A from 



S corr (k,E) = J d 3 r PA ( T )S™(k,E; 



p = p A (r)) , (31) 



where Pa( t ) is the nuclear density distribution and 
Scorrfa: E; p) is the correlation part of the spectral func- 
tion of uniform nuclear matter at density p. The corre- 
lation part of the nuclear matter spectral function can 
be easily singled out at zero-th order of CBF, being as- 
sociated to two hole-one particle intermediate states. At 
higher orders, however, one hole and two hole-one parti- 
cle states are coupled, and the identification of the cor- 
relation contributions becomes more involved. A full ac- 
count of the calcula tion of S^f.(k,E) can be found in 
(|Benhar et all Il994[) . 

The full LDA spectral function is written in the form 



S LDA (k,E) = S M f{KE) + S corr (k,E) 



(32) 



the spectroscopic factors Z n of Eq. ([3U|) being constrained 
by the normalization requirement 



d 3 k dE S LDA (k,E) = 1 



(33) 



A somewhat different implementatio n of LDA has also 
been proposed ( Van Neck et all Il995l ). Within this ap- 
proach the nuclear matter spectral function is only used 
at k > kp(r), kp(r) being the local Fermi momentum, 
whereas the correlation background at k < kp(r) is in- 
corporated in the generalized mean field contribution. 
Comparison between the resulting Oxygen momentum 
distribution and that obtained by Benhar et al. shows 
that they are in almost perfect agreement. 

The LDA scheme is based on the premise that short 
range nuclear dynamics are unaffected by surface and 



shell effects. The validity of this assumption is supported 
by the results of theoretical calculations of the nucleon 
momentum distribution 



i(k) = J dE [ZS p (h, E) + (A- Z)S n (k, E)] 



(34) 



showing that for A>4 the quantity n(k)/A becomes 
nearly independent of A at large |k| (> 300 MeV/c). This 
feature, illustrated in Fig. |H suggests that the correlation 
part of the spectral function also scales with the target 
mass number, so that S^ ^(\i, E) can be used to approx- 
imate S corr (k, E) at finite A. 
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FIG. 4 Calculated momentum distributi on per nucleon in H , 
4 He, 16 Q and uniform nuclear matter (jBenhar et al ], 1 19931 ; 
ISchiavilla et al\ , \ 19861 ). 



A direct measurement of the correlation component 
of the spectral function of 12 C, from the (e,e'p) cross 
section at missing momentum and energy up to ~ 800 
MeV/c and ~ 200 MeV, respectively, has been r ecent! 
carrie d out by the JLab E97-006 Collaboration ( Roh 
|2004[ ) . The data from the preliminary analysis appear to 
be consistent with the theoretical predictions based on 
LDA. 



D. Contribution of inelastic processes 

The approach described in the previous sections is not 
limited to quasielastic processes. The tensor defined in 
Eqs. (fTS)) and (fT9]) describes electromagnetic transitions 
of the struck nucleon to any hadronic final state. 

To take into account the possible production of 
hadrons other than protons and neutrons one has to re- 
place Wi and ui2 given by Eqs. (|23)l and (|24|) with the 
inelastic nucleon structure functions extracted from the 
analysis of e lectron-proton and electro n-deuteron scat- 
tering data ( Bodek and Ritchie) . Il98lf ). The resulting 
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IA cross section can be written as in Eq. ([6]). the two 
nuclear structure fu nctions W\ and W2 being given by 
(jBenhar et qlj . ll997ft 



Widql.w) = J d 3 k dE ^ZS p (k,E) f^- 



<(|q|,5)+- 



l«£(|q|,£) |kxq| ; 



(35) 



and 



W 2 (|q|,cj) = J d 3 k dE^ZS p (k 7 E) 




w[{\cilu;)— — -1 



EkU) — k • q 



(36) 



where the dots denote the neutron contributions. 
Eqs. (f35|) and (l36l) are obt ained using the prescription of 
Eq. (USD (|de Forestl . [l983h to preserv e gauge invariance . 
Note that the standard expression ( Atwood and Westl . 
Il973l ). widely used in studies of nuclear effects in deep in- 
clastic scattering, can be recovered from the above equa- 
tions replacing ui — > u and Ek — * Ma — En. 



0.005 




FIG. 5 Inclusive electron scattering cross section at E e — 
3.595GeV and 9 — 30°. The data poi nts represent the e xtrap- 
olated nuclear matter cross section (|Dav et all I1989T ) while 
the solid and dashed lines show the results of IA calculations 
carried out with and without inclusion of the inelastic contri- 
butions, respectively (jBenhar et aU[l99lft . 

As an example, Fig. [5] shows the quasi-elastic (dashed 
line) and total (solid line) inclusive cross sections of uni- 
form nuclear matter, at beam energy E e = 3.595 GeV 
and scattering angle 8 = 30°, evaluated using a phe- 
nomenol ogical fit of the nucleon st ructure functions 



and wg (jBodek and Ritchie!. 1 1981ft and the above men 



The data show that the transition from the quasi elas- 
tic to the inelastic regime, including resonant and nonres- 
onant pion production as well as deep inelastic processes, 
is a smooth one, thus suggesting the possibility of a uni- 
fied representation. 

The approach based on NMBT and the IA yields a 
good description of the measured cross section at energy 
loss uj > 1 GcV, corresponding to 1 S 1.3 (note that in 
the kinematics of Fig. [5] the top of the quasi free bump 
corresponds tow = loqe ~ 1.4 GeV). On the other hand, 
the data at lower energy loss are largely underestimated. 

The failure of IA calculations to explain the measured 
cross sections at u <c loqe has long been recognized, and 
confirmed by a number of theoretical studi es, carried out 
using highly realistic spect r al functions dBenhar et al 



1989[ ICiofi degli Atti et all Il992t iMeier-Haiduk et al. 
19831 ). see e.g. Fig. [SI It has to be ascribed to FSI be- 



tween the struck nucleon and the spectator particles, that 
move strength from the region of the quasi free bump to 
the low w tail. This mechanism will be analyzed in the 
next Section. 




200 220 240 260 

a (MeV) 

FIG. 6 Inclusive electron scattering cross section at E e = 
7.26GeV a n d = 8° for 3 He. The data points are from 
(|Dav et all 11979ft . the solid line show s the IA calculation 
based on the 3 He spectral function (Mcicr-H aiduk et all 
1 19831 ) . Approximate values for the scaling variable x are in- 
dicated on top. 



In conclusion, NMBT and the IA provide a consistent 
and computationally viable approach, yielding a quanti- 
tative description of the data in both the quasi elastic 
and inelastic regime, with the only exception of the re- 
gion of very low energy loss. Theoretical studies in which 
nuclear binding effects are included using realistic spec- 
tral functions also provide a quantitative account of the 
size and density depende nce of the European Muon Col- 
laboration (EMC) effect (|Benhar et q^J . Il997l fl999ft . 



tioned spectral function (jBenhar et all Il989ft 



E. DifFerent implementations of the IA scheme 

In spite of the fact that the basic assumptions un- 
derlying IA can be unambiguously stated, in the litera- 
ture one finds two different definitions of the IA inclusive 



9 



cross section: the one discussed in Sec. III. At involving 
the target spectral function, and another one, written in 
terms of the target mome ntum distribution of Eq. |34|) 

(jRinat and Taraginl . flggl . 

It has been shown ( Benhar et "oil. l200lh that the defi- 
nition in terms of the spectral function follows from min- 
imal use of the assumptions involved in the IA scheme, 
and correctly takes into account the correlation between 
momentum and removal energy of the participating con- 
stituent. On the other hand, a more extended use of the 
same assumptions leads to a definition of the IA cross 
section in which the nucleon spectral function is written 
in the approximated form 



S(k,E) = -n(k)6 



|k| 2 
2m 



> (37) 



so that the information on the target removal energy dis- 
tribution is lost. 

Figure [7] shows a comparison between the response 
function of a nonrelativistic model of uniform nuclear 
matter 3 obtained from the full spectral function (solid 
line) and the approximation of Eq. (|37p (dashed line). 



the average nucleon removal energy 



E = / d 3 kdE ES(k, E) = 62 MeV . 



(38) 



This feature is illustrated by the diamonds of Fig. 
that show the response obtained replacing |k| 2 /2m — > 
(|k| 2 /2m — E) in the argument of the energy conserving 
5- function of Eq. (f37j) . The results of this calculation 
turn out to be much closer to those derived in the previ- 
ous sections. 

In addition to the shift in the position of the peak, 
the dashed and solid curves sizably differ at low energy 
transfer, where the response obtained using the momen- 
tum distribution is much larger. Obviously, to identify 
corrections to the response arising from mechanisms not 
included in the IA picture, one has to start from the def- 
inition involving the minimal set of approximations. The 
results of Fig. [JJ show that a quantitative understand- 
ing of the effects of FSI, which are known to dominate 
the inclusive cross section at low requires the use of 
the spectral function in the calculation of the IA cross 
section. 



2.5 



> 

£ 1-5 
'a 

H l.o 

K 



0.5 



0.0 



_, , 1 , , , , 1 , , , 

; |q| = 5 fm-\/f 


■ > ■ < i ■ > ■ >j 

■>v Nuclear 
\ \ Matter 


i 0/ 




' 0/ 




— 1 ©/ 


.... 'fVvf^ 



0.2 



0.4 0.6 
a [GeV] 



1.0 



FIG. 7 Response of uniform nuclear matter at equilibrium 
density and |q|=5 fin -1 . The solid and das hed lines have 
been obtained using the full spectral function ijBenhar et all 
I1989T ) and the approximation of Eq. (|37[) . respectively. The 
diamonds represent the res ults obtained when shifting the 
dashed line by E = 62 MeV (|Benhar et all l200ll) . 

The two responses have similar shape, their width be- 
ing dictated by the momentum distribution. However, 
they are shifted with respect to one another. The peak of 
the dashed curve is located at energy loss lu ~ |q| 2 /2m, 
corresponding to elastic scattering off a free stationary 
nucleon, whereas the solid line, due to the removal energy 
distribution described by the spectral function, peaks at 
significantly larger energy. The shift is roughly given by 



3 The response function shown in Fig. [7] is proportional to the 
inclusive cross section in the case of scattering of a scalar probe. 
The generalization to the electromagnetic target tensor of Eq. ifjTJ 
is straightforward. 



III. FINAL STATE INTERACTIONS 

The existence of strong FSI in quasi-elastic scattering 
has long been experimentally established. The results of 
a number of (e, e'p) measurements, covering the kinemat- 
ical domain corresponding to 0.5 ^ Q 2 £ 8.0 (GeV/c) 2 
(lAbbott et all Il998l: iGarino et all Il992t iGarrow et all 



l2002HO'Neill et aZ.l . ll995l ; lRohe et adl2005T ). clearly show 
that the flux of outgoing protons is strongly suppressed 
with respect to the I A predictions. The observed attenu- 
ation ranges from 20-40 % in Carbon to 50-70 % in Gold. 

Being only sensitive to rescattering processes taking 
place within a distance ~ l/|q| of the electromagnetic 
vertex, the inclusive cross section at high momentum 
transfer in general is largely unaffected by FSI. However, 
the systematic discrepancies between data and the re- 
sults of highly accurate IA calculations indicate that the 
effects of FSI can become appreciable, indeed dominant, 
in the low to region, where the cross sections produced by 
IA calculations become very small. As IA cross sections 
in this region are most sensitive to the high momentum 
and high removal energy tails of the nuclear spectral func- 
tion 5 a quantitative understanding of FSI is required to 
unambiguously identify correlation effects. 



Note that for spectral functions describing correlated nucleons E 
is much larger than the average of single-particle energies. 
5 When ignoring the E-distribution of S(k, E) the cross section at 
low u) is related to the strength at large k. However, the spectral 
function exhibits a strong correlation between large |k| and large 
E. For example, in nuclear matter at equilibrium density, more 
than 50% of the strength at |k| = 1.5/m — 1 (just above the Fermi 
surface) is located at E > 80MeV. 
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In inclusive processes FSI has two effects: i) an energy 
shift of the cross section, due to the fact that the struck 
nucleon moves in the average potential generated by the 
spectator particles and ii) a redistribution of the strength, 
leading to the quenching of the quasielastic peak and an 
enhancement of the tails, as a consequence of NN scat- 
tering processes coupling the one particle-one hole final 
state to more complex n-particle n-hole configurations. 

Early attempts to include FSI effects were based on the 
Green's function formalism and multiple scattering the- 
ory, leading to a descriptio n of the dynamic s in ter ms of a 
complex optical potential ( Horikawa et aZ.I . Il980h . How- 
ever, while providing a computationally practical scheme 
to account for the loss of flux in the one-nucleon removal 
channel, the optical potential model employed relies on 
the mean field picture of the nucleus, and does not in- 
clude the effect of dynamical NN correlations. 

A similar approach has been adopted ( Chinn et all 
11989ft to analyze the longitudinal and transverse re- 
sponses of Eq. ([7]) and investigate the possible impor- 
tance of relativistic effects using Dirac bound-state wave 
functions and optical potentials. Although the authors 
suggest that relativity may play an important role in sup- 
pressing the response functions, particularly Rl, the in- 
terpretation of their results is hindered by the large un- 
certainty associated with relativistic descriptions of nu- 
clear dynamics and the appearance of the deep poten- 
tials, driven by the attempt to describe correlated sys- 
tems using an independent particle model. Systematic 
studies of relativistic effects carried out within NMBT 
show that they are indeed rather small. For example, 
when using the relativistic kinetic energy operator and 
when including boost corrections to the NN potential, a 
change of the binding energy of nuclear mat ter at equilib- 
rium density of ~10% only has been found ( Akmal et all 
11998ft . 

A different approach, based on NMBT and a gener- 
alizatio n of Glauber th eory of high energy proton scat- 
tering (iGlauberl . 1959ft was proposed in the early 90's 
( Benhar et all 1991 ). This treatment of FSI, generally 
referred to as the Correlated Glauber Approximation 
(CGA), rests on the premises that i) the struck nucleon 
moves along a straight trajectory with constant velocity 
(cikonal approximation), and ii) the spectator nucleons 
are seen by the struck particle as a collection of fixed 
scattering centers (frozen approximation). 

Under the assumptions given above the expectation 
value of the propagator of the struck nucleon in the tar- 
get ground state can be written in the factorized form 
(jPetraki et all . 12003ft 



with 



Ck+q(t) 



(39) 



where (t) is the free space propagator, while FSI 

effects are described by the quantity (R = (ri, . . . , r^)) 



U^(R-t) = I^c-^^./o^+.d^+w'l) _ (41) 



In Eq. (HTj) Vij = Yi—Yj and Fk+ q (|r|) is the coordinate- 
space t-matrix, related to the NN scattering amplitude 
at incident momentum k + q and momentum transfer p, 
Aq(p) by 



2?r 

r k+q (r) = 



"m / (^ PMk - (p) 



(42) 



At large |q|, k + q w q and the eikonal propagator 
of Eq. pij) becomes a function of t and the momentum 
transfer only. 

The scattering amplitude extracted from the measured 
NN cross sections is generally parameterized in the form 

A^p) = »M a NN (l - ia NN )e-^»^ 2 , (43) 

where o^n and oinn denote the total cross section and 
the ratio between the real and the imaginary part of the 
amplitude, respectively, while the slope parameter (3nn 
is related to the range of the interaction. 
The quantity 



^) = <o||cCW)llo) 



(44) 



measures the probability that the struck nucleon does 
not undergo rescattering process es during a t ime t after 
the electromagnetic interaction ( Rohe et all l2005h . In 
absence of FSI, i.e. for vanishing T q , Pq(t) = 1. 

Note that P(t) is trivially related to the nuclear trans- 
parency T n , m e asured in coincidence (e, e ' p) experiments 
(lAbbott et all 11998): iGarino et all Il992t iGarrow et all 



l2002t lO'Neill et all Il995l : iRohe et al.l . T2005ft . through 



T q = lim P q (t). 



(45) 



It is very important to realize that, as shown by 
Eqs. (|40p and (|4Tj) . the probability that a rescattering 
process occurs is not simply dictated by the nuclear den- 
sity distribution pA(jj), the probability of finding a spec- 
tator at position Yj. The rescattering probability de- 
pends upon the joint probability of finding the struck 
particle at position r j and & spectator at position Yj , that 
can be written in the form 



P (2) (ri,rj) = PA{Y l )pA{Yj)g{Y t ,Y J ) 



(46) 



Due to the strongly repulsive nature of nuclear inter- 
actions at short range, p( 2 '(ri,Tj) is dominated by NN 
correlations, whose behavior is described by the correla- 
tion function g{Yi,Yj). The results of numerical calcu- 
lations carried out within NMBT yield g(Yi,Yj) <C 1 at 
\Yij \ < 1 fm. This feature is illustrated in Fig. [SI showing 
the average p^ 2 - 1 (r^, rj), defined as 



\UtfZ(R;t)\0) 



(40) 



(Ty) = 2 J d 3 R l3 pW^rj) 



(47) 
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with Ry = (rj + Yj)/2, evaluated for both NN and pp 
pairs with and without inclusion of dynamical correlation 
effects. 




FIG. 8 Averaged distribution functions, defined by Eq. (|47|l. 
of NN (solid lines) and pp (dashed lines) pairs in 4 He, ie O 
and uniform nuclear matter at equilibrium density (NM). The 
corresponding quantities obtained neglecting correlation ef- 
fects are shown by th e dot-dash and dotted lines, respectively 
(jBenhar et gill 1993t h 

From Eq. (gT]) it follows that within the CGA the en- 
ergy shift and the redistribution of the inclusive strength 
are driven by the real and the imaginary part of the NN 
scattering amplitude, respectively. However, at large q 
the imaginary part of T q , which gives rise to the real part 
of U^ SI , is dominant. 

Neglecting the contribution of the real part of Aq alto- 
gether, the CGA inclusive cross section can be written as 
a convolution integral, involving the cross section evalu- 
ated within the IA, i.e. using Eqs. ^ and ([TO]) , and 
a folding function embodying FSI effects: 



da 



dfl„i du> 



duj' 



da 



dQp/duj' J IA J 



/ q (w - w') , (48) 



/ q (w) being defined as 



dt 
2^ 



/ q (w) = <5^)v/T q + 

= i(u)^ + f q (w) , 
and normalized according to 



+oo 



dufqibj) = 1 



(49) 



(50) 



The proceeding equations show that the strength of FSI 
is governed by both T q and the width of F q (w). In ab- 
sence of FSI Ug (t) = 1, implying in turn T q = 1 and 
/ q (w) -» 5(oj). 

In principle, the real part of the NN scattering ampli- 
tude can be explicitly included in Eq. (|41j) and treated 



on the same footing as the imaginary part. However, its 
effect turns out to be appreciable only at t ~ 0, when the 
attenuation produced by the imaginary part is weak. The 
results of numerical calculations show that an approxi- 
mate treatment, based on the use of a time independent 
optical potential, is adequate to de scribe the energy shif t 
produced by the real part of r q ( Benhar et aLl . Tl994h . 
whose size of ~10 MeV is to be compared to a typical 
electron energy loss of a few hundred MeV. 

The shape of the folding function is mainly determined 
by the total NN cross section &nn- In the energy re- 
gion relevant to scattering of few GeV electrons, er/vjv 
is dominated by the contribution of inelastic processes 
and nearly independent of energy. As a consequence, the 
scattering amplitude of Eq. (|4"3"]) grows linearly with q 
and both Uq SI (t) and the folding function / q become 
independent of q. 
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FIG. 9 Folding functions _F q of Eq. (|49|) . calculated in uni- 
form nuclear matter at equilibrium density. The solid line 
and the diamonds correspond t o momentum transfer |q| = 
2.2 and 3.4 GeV/c, respectively (|Benhaii Il999h . 



This feature is illustrated in Fig. [5] showing that the 
Fq calculated in uniform nuclear matter at equilibrium 
density at momentum transfers |q| = 2.2 and 3.4 GeV/c 
are nearly identical. 

Dynamical NN correlations also affect the shape of the 
folding function. Inclusion of correlations through the 
distribution function g(ji,rj) results in a strong quench- 
ing of the tails and an enhancement of the peak of Fq, 
leading to a significant suppression of FSI effects. 

The effect of FSI is illustrated in Fig. [TO] showing the 
inclusive cross section of uniform nuclear matter at a 
beam energy E e = 3.595 GeV and a scattering angle 
9 = 30°, corresponding to momentum transfer q| ~ 2 
GeV/c. 

Comparison between theory and the data in Fig. [TO] 
clearly shows that at w < 1.1 GeV, where quasiclastic 
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scattering dominates 6 and which correspond to x > 1, 
x = Q 2 /2muj being the Bjorken scaling variable, FSI 
effects are large and must be taken into account. The 
results obtained within the CGA are in good agreement 
with the data in the region co > 800 MeV, i.e. for x < 1.8, 
while at higher x the experimental cross section is largely 
overestimated. The dashed line has been obtained ne- 
glecting the effect of dynamical correlations on the distri- 
bution function g(ri,Tj). Comparison between the solid 
and dashed lines provides a measure of the quenching of 
FSI due to NN correlations. 
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FIG. 10 Inclusive electron scattering cross section at E e = 
3.595GeV and 9 — 30°. The data poi nts represent the extrap- 
olated nuclear matter cross section (|Dav et all 1 19891 ). while 
the solid and dashed lines show the results obtained in CGA 
including FSI effects, with and without taking into account 
correlation effects. For compari son, the IA c r oss se ction is 
also shown by the dot-dash line (|Benhar et all Il99ll ) . 



The ability of the CGA to provide a quantitative un- 
derstanding of FSI in the region x < 2 is further illus- 
trated in Fig. QT] showing the cross section ratio 



I? 



da(e + 56 Fe 



X) 2 



da(e + 2 H 



X) 56 



(51) 



at E e = 3.595GeV and 9 = 25°. Note that R of Eq. d5TJ) 
is only defined up to y ~-700 MeV/c, corresponding to 
x = 2, the kinematical limit for inclusive scattering off 
an A=2 target (for the definition of y see Sec. IVII|) . 

The solid line in Fig. [IT] corresponds to the full CGA 
calculation, providing a good description of the experi- 
ments over the whole range of y, whereas the IA results, 
represented by the dashed line, lie well below the data 
at y < -200 MeV/c (x > 1.5). For comparison, Fig. [TT1 
also shows the results obtained using the approximate 



In the kinematics of Fig. 1101 inelastic processes only contribute 
~ 5% of the inclusive cross section at lu = 1.1 GeV, and become 
negligibly small at lower lo. 
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FIG. 11 Ratios of inclusive cross sections of iron and deu- 
terium at E e = 3.595GeV and 9 = 25°. Solid line: full cal- 
culation; dashed line: IA calculation, neglecting FSI in both 
iron and deuteron; dot-dash line: calculation ca rried out usin 
the ap proximate spectral function of Eq. (|37[1 (jBenhar et al 
Il995bh . 



spectral function of Eq. (|37|) . which turn out to largely 
overestimate the data at negative y. 

The results of Figs. [11] and [S] clearly rule out the in- 
terpretation of the behavior of the cross section ratio at 
y ^ —200 MeV/c in te rms of scatter i ng off strongly cor- 
related nucleon pairs (|Egivan et all I2003D . This inter- 
pretation in fact assumes the validity of the IA picture, 
which is known to fail at large negative y, and does not 
take into account the large effect of FSI. 

Notwithstanding its success in describing the existing 
inclusive data at large negative y, the CGA appears to 
consistently overestimate FSI effects at larger —y. As the 
validity of the eikonal approximation is well established 
in the kinematical region apposite to scattering of few 
GeV electrons, possible corrections to the CGA scheme 
are likely to be ascribable either to modifications of the 
NN scattering amplitude or to the inadequacy of the ap- 
proximations leading to the convolution expression for 
the cross section. 

It has been pointed out ( Benhar et aZ.I . [l99ll ) that the 
use of the free-space amplitude to describe NN scattering 
in the nuclear medium may be questionable. Pauli block- 
ing and dispersive cor rections are known to be importan t 
at moderate energies ( Pandharipande and Pieper[|l992t h 
However, their effects on the calculated inclusive cross 
section have been found to be small in the kinematical 
region corres ponding to |q| ^ 2 Ge V/c, and decrease as 
q increases (jBenhar et aUll995a| ). Corrections to the 
amplitude associated with its extrapolation to off-shell 
energies are also expecte d to be small at |q| > 2 GeV/c 
(|Benhar and Liutil . Il996h . 

Modifications of the free-space NN cross section may 
also originate from th e internal structure of the nuc leon. 
It has been suggested (iBrodskvi . 119821 : iMuellerl . Il982h that 
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clastic scattering on a nuclcon at high momentum trans- 
fer can only occur if the nucleon is found in the Fock 
state having the fewest number of constituents, so that 
the momentum transfer can be most effectively shared 
among them. This state is very compact, its size being 
proportional to l/|q|, and therefore interacts weakly with 
the nuclear medium. Within this picture a nuclcon, af- 
ter absorbing a large momentum, travels through nuclear 
matter experiencing very little FSI, i.e. exhibits color 
transparency (CT), before it evolves back to its standard 
configuration with a characteristic timescale. 

CT may be particularly relevant to the analysis of in- 
clusive electron- nucleus scattering at x > 1, where elastic 
scattering is the dominant reaction mechanism, since it 
leads to a significant quenching of FSI. In fact, the influ- 
ence of CT is expected to be much larger for (e, e') than 
for (e, e'p); in (e, e') FSI occur mainly very close to the 
electromagnetic vertex, at a distance of less than ~ 1/ |q| , 
where the compact configuration has not yet evolved back 
to the ordinary proton. It has been suggested that the 
modification of the NN scattering amplitude due to on- 
set of CT may explain the failure o f CCA to reproduce 
the data in the region of very low a ; (IBenhar et aZ.I . [l99lL 
11994 iBenhar and Pandharipandd . [l993lh 

However, the recent measurements show no enhance- 
ment of the nuclear t ransparency u p to Q 2 ~ 8 (GeV/c) 2 
(jDutta et all 120031 ; iGarrow et all . |2002( ) and seem to 
rule out observable CT effects in (e,e'p) at beam en- 
ergies of few GeV. Therefore, the excellent agreement 
between the measured inclusive cr oss sections a n d the 
resul t s of theoretical calculation (IBenhar et all Il99l[ 
Il994t IBenhar and Pandhariparidel Il993[ ) carried out us- 
i ng CGA and the quantum diffusion model of CT 
( Farrar et all fl988h may be accidental. 

An i mproved version of CGA has been recently pro- 
posed ( Petraki et all 120031 ). Within this approach the 
initial momentum of the struck nucleon, which is aver- 
aged over in CGA, is explicitly taken into account. As 
a result, one goes beyond the simple convolution form 
of the inclusive cross section and gets a generalized fold- 
ing function, depending on both the momentum transfer 
q and the initial momentum k. Numerical calculations 
of the nonrelativistic response of uniform nuclear mat- 
ter at 1 ;$ |q| ;$ 2 GeV/c show that the inclusion of 
this additional momentum dependence leads to a sizable 
quenching of the low energy loss tail of the inclusive cross 
section, with respect to the predictions of CGA. 

Even though they cover a limited kinematical range 
and have been obtained using a somewhat oversimplified 
model, the available results suggest that a better treat- 
ment of the momentum distribution of the struck nucleon 
may improve the agreement between theory and data in 
the region of x > 2, where CGA begins to fail. A sys- 
tematic study of the dependence of the rescattering prob- 
ability on the initial momentum, based on a relativisti- 
cally consistent formalis m, is presently being carried out 
(IBenhar and Sickl . l2006lh 

A different approach to describe FSI in inclu- 



sive electron-nucleus scattering, based on the rela- 
tivistic generalization of the Gcrsch-Rodrigucz-Smith 
(GRS) l/|q| e xpansion of the r esponse of many- 
fermion systems dGersch et all [l973h . has been proposed 
(jGurvitz and Rinatl 12002 ). Numerical studies of the 



cross section of 4 He ( Viviani et al I l2003h . carried out 
using realistic wave functions to compute the ( 4 He\ 3 He') 
and l^He^H^ tra nsition matrix elements and a mo del 
spectral function ( Ciofi degli Atti and Simula! Il994h to 
describe the many-body breakup channels, show that the 
generalized GRS expansion pr ovides a fai r ly goo d overall 
description of the SLAC data ( Dav et oZl . ll993h at beam 
energy of 3.6 GeV. 



IV. EXPERIMENTS 

Studies of inclusive electron-nucleus scattering have 
been performed at a number of facilities. The bulk of 
the currently available data (at momentum transfers |q| 
on the order of 500 McV/c and above upon which we fo- 
cus here) has been generated at Saclay (France), Bates- 
MIT (Boston), Jefferson Lab (Newport News) and SLAC 
(Stanford). Isolated data sets are available from other fa- 
cilities (see next section). 

Without exception, magnetic spectrometers have been 
used to momentum-analyze the scattered electrons. The 
magnet arrangements are varied, and at the four facilities 
mentioned above were D(ipole), DD, Q(uadrupole)QQD 
and QQDDQ systems, respectively. These spectrome- 
ters had a maximum bending momentum of 0.6, 0.9, 7 
and 8 GeV/c , respectively and solid angles between 4 



and 6.8 msr (IBertozzi et all Il979t iLeconte et all . Il980t 



iMo and Pecld . 1 1 96.4 1 Yanl . 1 1 9 9 7t ) . As an example, we show 
in Figs. 1121131 the JLab High Momentum Spectrometer 
(HMS) of Hall C and its associated focal-plane detector 
setup. 




FIG. 12 HMS spectrometer in hall C at Jefferson Lab with 4 
magnets arranged as QQQD and the heavily shielded detector 
hut. 



In most cases the electron track in the focal plane was 
reconstructed using 2 to 4 planes of multiwire chambers, 
and the trigger involved the use of fast scintillator detcc- 
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tors, often 2 planes of segmented paddles. Typically the 
focal plane detector also included a Cerenkov detector 
needed to separate electrons from slower particles such 
as pions. At the high-energy facilities a total absorp- 
tion counter provided an additional, cruder, energy mea- 
surement to assist in the discrimination of electrons from 
other, often more numerous, charged particles. 



due to detector acceptance and ine fficiencies and the true 
^-dep endence of the cross section ( Crannell and Suelzld . 
1196 



drift 
chambers 



o 
o" 



Cerenkov 



scintillators 



2.4 m 



lead glas: 



FIG. 13 HMS detector package, consisting of multiple 
planes of multiwire drift chambers, fast scintillators, a sub- 
atmospheric gas Cerenkov counter and a lead glass total ab- 
sorption counter. 



The data acquisition systems of the earlier, lower- 
energy facilities were rather restricted because of com- 
puter speed, memory and high data rates. Consequently 
much of the data was pre-processed online and stored in 
spectra; this allowed little off-line analysis of the data. 
The data acquisition at the higher energy facilities in- 
cluded event-mode acquisition where all the coordinates, 
times and amplitudes are registered for further off-line 
analysis. Event-mode acquisition has important advan- 
tages in that during the course of the off-line analysis 
shortcomings of the acquisition hardware can be discov- 
ered and (often) corrected. This typically leads to much 
more reliable cross sections. 

In many cases the product of effective solid angle and 
detector efficiency is not easily established. While the 
efficiency of the detector system can be determined by 
exploiting the redundancy of the detector elements, the 
effective solid angle of spectrometers (not using a solid- 
angle defining slit of well-known geometry) is difficult 
to establish. For that reason many of the experiments, 
in particular the ones at the higher-energy facilities, use 
elastic scattering from a liquid hydrogen target as check; 
the e-p cross sections are well known. The data on hy- 
drogen often are also used to determine offsets in the 
scattered electron momenta and angles as reconstructed 
from positions and angles of the electrons in the detec- 
tors. 

The inclusive quasielastic spectrum extends over a mo- 
mentum range that exceeds (in most cases) the momen- 
tum acceptance of the spectrometers. Consequently data 
taking entails a sequence of runs at stepped values of the 
central spectrometer momentum. The resulting set of 
spectra, if taken with sufficient overlap, can then be used 
to determine both the dependence in Sp = Ap/p centr ai 



An experimental aspect peculiar to inclusive scattering 
is the need to accurately know the dispersion of the spec- 
trometer. As the cross section is ultimately determined 
for bins in energy loss, the accuracy of the differential mo- 
mentum dispersion enters directly in the accuracy of the 
cross sections. The dispersion can be determined using 
focal-plane scans with an elastic peak, provided that the 
relation between the integral over Bdl and the quantity 
measured — typically the B-field at one point — is linear 
and has no offset at zero momentum. For spectrometers 
with high resolution the dispersion can be determined us- 
ing elastic scattering off nuclei with different mass. This 
determination of a precise dispersion often has not been 
given the appropriate attention. 

After correcting the data for the various shortcomings 
of the spectrometers and detector packages (energy- and 
angle-offsets, detector inefficiencies, ...) one major cor- 
rection is the removal of radiative effects from external 
and internal Brcmsstrahlung. Whi le the theory for ra - 
diativc effects is well under control ( Mo and Tsai 119691) . 
the specifics of quasi-elastic scattering lead to the need 
for extensive modeling. In quasi-elastic scattering, the 
energy loss can be very large as can be the energy loss 
due to radiative effects. Since radiation can be emit- 
ted before or after the quasi-elastic scattering process, 
the calculation of radiative effects requires knowledge of 
the quasi-elastic structure function for all incident elec- 
tron beam energies down to E e — u) max , where E e is the 
incident energy and uj max is the largest energy loss of in- 
terest. Collecting the required quasi-elastic data at the 
corresponding lower incident electron energies is often 
too time consuming. As a consequence, only few such 
data or none at all are taken. In order to obtain the 
quasi-elastic cross section at all energies extensive mod- 
cling of the data in the region covered by experiment is 
needed, such that a reliable extrapolation to the other 
regions is feasible. With the discovery of y-scaling (see 
Sec lVIip this modeling has been greatly facilitated, as the 
scaling function F(y) represents a quantity that is only 
very weakly dependent on energy and angle and can be 
interpolated over a large range. This unfortunately does 
not apply to the A-region, which at large uj overlaps with 
the quasi-elastic peak. Modeling of the A contribution 
can be faci litated when using t he approach based on su- 
perscaling ( Amaro et all |2005| ) . 

Despite the difficulties introduced by the radiative cor- 
rections, the impact on the overall accuracy is relatively 
small. For cross sections that do not depend steeply on u 
the effect of the radiative corrections is small; the uncer- 
tainty introduced is also small if the model cross section 
is not grossly deficient. 

There is another correction — not always considered 
— that is unique to cross section measurements of con- 
tinuum spectra at large energy losses. One must account 
for the contributions of electrons with momenta outside 
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the momentum acceptance of the spectrometer which in- 
teract with part of the spectrometer yoke or vacuum 
chamber, and are rescattered into the detector. Usu- 
ally the trajectory of these electrons cannot be traced 
back through the spectrometer as their momentum is not 
known. This process can lead to background, which in 
general is in the %-region, but can become problematic 
when e.g. an L/T-separation is performed. In this case 
the effect is greatly enhanced in magnitude if the L- or 
T-contribution is small as compared to the other one, 
as happens at the larger momentum transfers where the 
longitudinal response gives a small contribution to the 
cross section. Such rescattering events can become es- 
pecially important if the experiment uses long (liquid or 
gas) targets that extend beyond the acceptance of the 
spectrometer; in this case electrons can enter the spec- 
trometer through the solid-angle defining collimator and 
scatter from the vacuum chamber or pole pieces. 

Experimentally, provisions against scattering from the 
yoke (but not the pole pieces) include the use of baffles 
on the high/low mom entum wall of the vacuum chamber 
( Leconte et all Il980l ) . In some cases these backgrounds 
also have been subtracted by Monte-Carlo simula tions of 
the p aths of electrons through the spectrometer ( Danell . 
1990). They have been measured using a solid target dis- 
placed along the beam direct i on, an d corrections up to 
10% have been found (jDanell . [199(1 . The most reliable 
method to remove these contributions has been accom- 
plished via a secondary measurement of the energy of the 
electron in the fo c al pla ne, using some total-absorption 
counter ( Jourdanl fl996). Once one knows with reason- 
able accuracy the energy of the detected particles, they 
can be traced back and one can find out whether they 
come directly from the target or not. For the QQDDQ 
8GeV spectrometer at SLAC an acceptance-reducing col- 
limator placed in between the two dipoles allowed an ex- 
perimental determination of the fractional contribution 
of rescattered events. At large energy loss this f ractio n 
has been found to be quite substantial ( Jourdanl . Il996l ). 



On the whole, this correction for rescattering has 
been given too little attention, particularly at the lower 
energies where their contribution is more important 
and where secondary energy-measurements using total- 
absorption counters are not very effective. 

One further correction that has to be applied to the 
data is the subtraction of electrons resulting from e + e~ 
pair production in the target. The e~ from this charge- 
symmetric process give a significant contribution to the 
inclusive cross section, particularly at large scattering an- 
gles and small energy loss to where the quasi-elastic cross 
section is very small. The effect of this contribution 
is usually handled by subtracting the spectrum, mea- 
sured under identical kinematic conditions, of the pair 
produced e + . 



V. DATA 

The modern era of experiments dedicated to quasielas- 
tic electron-nucleus scattering be gan in 1969 wit h a se- 
ries of data sets from Kharkov (IDementii et al. I. ll!Mi!t 
iTitov and Stepulal . Il972l : iTitov et aZl . ll969l . ll97ll . Il974h : 
a first systematic survey of quasielastic electron scatter- 
ing from a range of nuclei (He to Pb) ( McCarthy et all 



19741 ) was made 



119761 ; iMoniz et all Il97ll ; IWhitnev etal. 
at Stanford. The first attempt to separate the re- 
sponse functions was performed at Stanford when data 
for 40 C a and 48 Ca were collected at constant values 
of the three-momentum trans fer ( Zimmerman! . Il969t 
IZimmerman and Yearianl Il976l ). 

At about the same time quasielastic electron scattering 
from 12 C was being studied at the Harvard CEA at ener- 
gies ranging from 1 to 4 GeV and angles from 8.5° to 18° 
using a half-quadrupole spectrometer ( Stanfield et all 
At DESY inclusive electron scattering 



119711 



data 



in the quasiela stic re gion w as studied ( Heimlich! Il973t 
iHeimlich et~al\ . Il974t IZeilerl . Il973l ) from 6 Li and 12 C at 

energies up to 2.7 GeV and angles as large as 15 degrees. 

The high energy beam at SLAC was exploited in 
the mid-1970's to measure quasielastic cross sections at 
very high momen tum transfer (~ 4 (GeV/c) 2 ) on 3 He 



( Day et all . I1979T ) . A considerable body of quasielas- 
tic data was also measured at SLAC in the mid- 



1980's (iBaranl. Il989t iBaran et all 19881 


Chen. 199C; 


Day et al. 19931. 1987; Meziani et all 1992 


; Pottcrvcld, 


19891; ISealock et all Il989h using the high-intensity, low 



energy (<5GeV) beam from the NPAS injector. 

Starting in the late 1970's a series experiments 
began to measure data for a wide range of nuclei ( 2 H 
to U), often as part of a program to separate the 



structure functions, both at Bates (lAltemus et al. 


1980 


; iBlatchlev et all 1986; Deadv et al. 19861 


1983 


Dow 


, 19871 Dow et al. 1988; Dvtman et all 


1988 


Hotta et all 1984; O'Connell et all 1987; Quinn et al. 


1988; VanReden et al. 


. 19901; Yates et al. 1993 


Zimmerman et al.. 1978 


) and at Saclav (Barreau et al. 


1983; Gueve et al. 19991 iMarchand et all 


1985 


Meziani et all 19851 1984 Mougev et all 


1978 


Zghiche et aUll994) 



In 1996 an experiment at JLab, intended to extend the 
range in four-momentum transfer Q 2 and Bjorken scaling 
variable x, pro duced data in the quasielastic reg i on ou t 
to 5 (GeV/c) 2 (|Arringtonl . Il998t lArrington et aZ.I . Il999h . 



In the last year a new experiment using the full energy of 
Jlab (6 GeV) extended these measurements even further 
using 2 H , 3 He, 4 He, Be, C, Cu and Au targets. 

As most of the data is not published in numeri- 
cal form, we have prepared, as part of this review, a 
website Quasielastic Electron Nucleus Scattering Archive 
|http : / /faculty . V irginia . edu/ qes- archive. In this 
archive we have placed all radiatively unfolded cross sec- 
tion we could locate in tabular form. The site also gives 
the references and some details on the experiments. 

For some nuclei, quite a comprehensive set of cross sec- 
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FIG. 14 Kinematic range covered by the (e,e') data for 12 C, 
plotted as function of |q| and y (for definition see Sec. IVIIfl . 
The symbols encode the range of e covered by the data (bins 
of 0.2 in e). 



tions is available. Fig. [14] shows the range in |q| and y 
covered by the data for 12 C. The symbols incorporate 
information on e = l./(l. + 2q 2 /Q 2 tan 2 8/ 2), indicat- 
ing that only for the lower momentum transfers a large 
enough e-range is covered for reliable L/T-separations. 



VI. EXTRACTION AND USE OF NUCLEON FORM 
FACTORS 

We began our discussion of the quasielastic scattering 
process by describing it as the incoherent scattering of 
the high energy electrons from the bound and quasi-free 
nucleons. Any quantitative description of this process 
then depends on an accurate description of the elemen- 
tary electron-nucleon (proton and neutron) cross section. 

In this section we will briefly review the status of the 
form factors, with particular attention to the most re- 
cent developments. A comprehensive review is not called 
for a n d the interested reader shou l d con sult Refs. ( Gaol . 
120031: iHvde- Wright and de Jagerl. I2004D or the r ecent 
conference proceedings ( Davl . 20051 ; de Jagerl . l2004h . 

The form factor data, until recently, have been ob- 
tained, in the case of the proton, through a Rosenbluth 
separation which exploits of the linear dependence of the 
cross section on the polarization e of the virtual photon. 
In terms of exchange of a single virtual photon the cross 
section depends on the Sachs form factors Ge and Gm 



da 



G E 



1 



^i + 2rG 2 / tan 2 (0/2) 



(52) 



with r = Q 2 /AM 2 and ctns = <^MottE'/EQ. Rearranging 



this expression with e 1 = 1 
da e(l 



or 



dfl 



CTNS 



2(1 + r)tan 2 (6>/2) yields 
T G 2 M (Q 2 ) + eG 2 E (Q 2 ). (53) 



By making measurements at a fixed Q 2 and variable 
e(9 7 E ), the reduced cross section ctr can be fit with a 
straight line with slope G 2 E and intercept tG 2 m . The 
Rosenbluth technique, because of the dominance of 
G P M over G P E at large Q 2 , demands strict control over the 
kinematics, acceptances and radiative corrections over a 
large range of both incident and final electron energies 
and angles. 

The analysis of quasielastic scattering has utilized 
widely accepted parameterizations or models of the nu- 
cleon form factors, based on data collected, starting in 
the 1950's, through the 1990's. The data has been 
described by simple fits as well as sophisticated mod- 
els based on the Vector Dominance Model VDM 
persion relations and quark models, see dGao 
IHvde- Wright and de Jagerl |2004 iPerdrisat et al 
for a review. 

The lack of a free neutron target with which to 
study the neutron properties forces the experimental 
effort to turn to investigate the neutron bound in a 
nucleus. The deuteron, weakly bound and with a 
ground state that, in principle, is calculable, has been 
for unpolarized studies the nucleus of choice. Un- 
til the early 1990's the extraction of G E was done 




most successfully through either small angle elastic 
electron scattering from the deuteron (|Benaksas et al.. 


1964; Bumiller et all 1970t Drickev and Hand. 1962 


Galster et o/lll97ll Platchkov et all 19901; 


Simon et at. 


1981 


) or by quasielastic e-D scattering ( 


Bartel et at. 


1973 


;lHughes et all 1 19661 Il965l iLung et al. 


. 1993h. 



In the IA the elastic electron-deuteron cross section is 
the sum of proton and neutron responses with deuteron 
wave function weighting. The coherent nature of elastic 
scattering gives rise to an interference term between the 
neutron and proton response which allows the smaller 
G E contribution to be extracted. Still, the large proton 
contribution must be removed. Experiments have been 
able to achieve small statistical errors but remain very 
sensitive to deuteron wave function model leaving a sig- 
nificant residual dependence on the NN potential. 

Quasielastic e-D scattering provides a complementary 
approach to the extraction of G 7 ^. In th e IA model pi- 
oneer ed by Durand and McGee ( Durandl . Il96ll ; iMcGed . 
I1967D the cross section is incoherent sum of p and n cross 
section The extraction of G E requires both a Rosenbluth 
separation and the subtraction of the sizeable proton con- 
tribution. It suffers, unfortunately, from unfavorable er- 
ror propagation and a sensitivity to the deuteron struc- 
ture. 

G^j has been extracted from both inclusive and exclu- 
sive quasielastic scattering from the deuteron. G 7 ^ when 
determined from the ratio, 2 H(e, e'p)/ 2 H(e, e'n), always 
working near the top of the quasielastic peak, has the 
smallest uncertainties arising from FSI, MEC and details 
of the ground state wavefunction. For both the ratio 
technique and 2 H(e, e'njp the absolute efficiency of the 
neutron detector must be determined. Th e most precise 
data arc from MAMI (jKubon et all l2002h . recently G 7 ^ 
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was measured out to 5 (GeV/c) 2 dBrooks and Lachnietl . 
120051: lLachnietl . l2006h at JLab. 

It has been known for many years that the " small" nu- 
cleon electromagnetic form factors, G E and G P E at large 
Q 2 could be measured through spin-dependent elastic 
scattering from the nucleon, accomplished either through 
a measurement of the scattering asymme try of polarized 
electrons from a polarized nucleon target (IDombeyl . ll969t 
iDonnellv and Raskinl . 1 19861 : [Raskin and DonnelM Il989h 
or equivalently by measuring the polarization transferred 
to th e nucleon ( Akhiezer and Rekalol . fl974ll Arnold et all 
Il98lh . 

In the scattering of polarized electrons from a polarized 
target, an asymmetry appears in the elastic scattering 
cross section when the beam helicity is reversed due to 
the presence of a polarized piece, A, in addition to the 
unpolarizcd piece, S. The total cross section is 



a(h) = £ + hA; h = ±pb eam - 



The asymmetry 



0V 



A 



(54) 



(55) 



can be written schematically (a, b, c, and d are known 
kinematic factors) as 



.4 



a cos 9* (G M ) + b sin 9* cos <f>*G E G 



A I 



c(Gm) 2 +d(G E ) 2 



(56) 



where 9* and <E>* fix the target polarization axis. With 
the target polarization axis in the scattering plane and 
perpendicular to q, (9*,$* = 90°, 0°) the asymme- 
try Atl is proportional to GeGm- With the polar- 
ization axis in the scattering plane and parallel to q 
(9*,$* = 0°,0°), one measures the transverse asymme- 
try At, which in the case of a free nucleon would be 
insensitive to G M (depend simply on the kinematic fac- 
tors). In 3 Hc(e, e')X scattering the denominator contains 
contributions arising from both the protons and the ne u- 
trons and at modest Q 2 is sensitive to G M ( Gaol . 120031 ) . 

In elastic scattering of polarized electrons from a nu- 
cleon, the nucleon obtains (is transferred) a polariza- 
tion whose components, Pi (along the direction of the 
nucleon momentum) and Pt (perpendicular to the nu- 
cleon momentum) are proportional to G 2 M and GeGm 
respectively. Polarimeters are sensitive only to the per- 
pendicular polarization components so precession of the 
nucleon spin before the polarimeter in the magnetic field 
of the spectrometer (for the proton) or a dipolc (in- 
serted in the path of neutron) allows a measurement of 
the ratio Pt/Pi and the form factor ratio: Ge/Gm = 
— jf - tan(fl/2). The recoil polarization technique 

has allowe d precision measurements of Gg to nearly 6 
(GeV/c) 2 (iGavou et all 1200 ll I2002H Jones et aZ.l . l2000D . 

These data on the proton have generated a great deal 
of activity as they have revealed a major discrepancy 



between the Roscnbluth data and the polarization trans- 
fer data for G P E at Q 2 > 2GeV 2 /c 2 . A detailed dis- 
cussion of these topics is beyond the scope of this re- 
view; it is discussed in the recent papers dGa^l 120031: 
IHvde- Wright and de Jagerll2004 fPerdrisat et al[|2006[ ) 



Extraction of the neutron form factors using polariza- 
tion observables is complicated by the need to account for 
Fermi motion, MEC, and FSI, complications that are ab- 
sent when scattering from a proton target. Fortunately it 
has been found for the deuteron that in kinematics that 
emphasize quasi -free neutron knoc kout both the transfer 
polarization Pt ( Arenhoevell. 19871) and th e beam-target 
asymmetry A^P (jArenhoevel et all Il988l ) arc especially 
sensitive to G E and relatively insensitive to the NN po- 
tential describing the ground state of the deuteron and 
other reaction details. See Figure [TSl The utility of 
3 He as a polarized neutron arises from the fact that in 
the ground state (dominantly a spatially symmetric S 
wave) the proton spins cancel and the 3 He s pin is carried 
by the unpaired neutron. Calculations (jGolak et all 
120011: llshikawa et al. I. I1998D of the beam-target asymme- 
try from a polarized 3 He target have shown it to have 
only modest model dependence. 

Extraction of the neutron form factors from beam- 
target asymmetry measurements in inclusive scattering 
are hindered by the dominant contribution of the pro- 
ton, even with a polarized 3 He where the protons are 
responsible for just 10% of the 3 He polarization. This 
is especially true in th e case of the neutron charge 
form factor ( Jones et al I 119931 [j ones- Woodward et all . 
Il99lh though methods to minimize the role of the pro - 
ton have been proposed ( Ciofi degli Atti et all Il995h . 



G M has been extracted from inclusive polarized electron- 
polarized 3 He scatteri ng at Bate s and J efferson Lab out 
to_Q 2 =0.6(GcV/c) 2 dGao et all 11994 IXu et all |2000L 
120031) . Inclusive as ymmetry measurem ents on 3 He in 
the threshold region (jXiong et al 1 l200lh have been used 
to successfully test nonrelativistic Faddeev calculations 
which include MEC and FSI. 

Fortunately the development of polarized beams and 
targets have been able to leverage the utility of exclu- 
sive quasielastic scattering from bound neutrons in both 
3 He and 2 H providing precision data on the electric form 
factors of the neutron. Coincidence measurements allow 
one to avoid completely the subtraction of the dominant 
proton. 

Since the e arly wo rk on the n eutron at Bates 
(|Eden et all 11994 iJones et all Il993l ) through 
3 Hc(e, e')X and D(e, e'n)p respectively, the further 
development of high polarization beams and targets, 
together with high duty factor accelera tors, has im- 
prove d the data set (a n d outl oo k) for G^ dBecker et al 
19991; [Herberg et all Il999l; iMeverhoff et all 11994 : 



Ostrick et all KL999I; iPasschier et all Il999l; iRohe et al , 
19991 : IWarren et all l2004 IZhu et all l200ll ) through 
either D(e,e'n)p, 3 He(e, e'n) or 2 H(e, e'n)p and G\ t 
through 3 He(e, e')X dGao et all 11994 IXu et all |2000L 
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FIG. 16 Re sult from the global an alysis of the L-T (Rosen- 
bluth) data (|Arrington et all 120061 ') . Open (full) points are 
obtained before (after) two-photon exchange corrections. The 
dashed line is a fit to the polarization-transfer data. 



FIG. 15 A calculation (|Arenhoevel et all 1 19881 . 1 19891 . Il992f ) 

of the electron-deuteron vector asymmetry A^ d at Q 2 = 0.5 
(GeV/c) 2 which demonstrates the insensitivity of A^ d to the 
reaction model and subnuclear degrees of freedom at angles 
near 180 deg. 
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to understand the discrepancy between 
transfer and R osenbluth resu lts for 
inve stigators ( Blunden et all 



Chen et al . 20041; Guichon and Vanderhaeghenl . 



2003; 



2003; 



Rekalo and Tomasi-Gustafssonl 2004 ) have explored 



the possibility of two-photon exchange corrections. 
While only incomple te calculations exist , the results of 
Ref. (|Blunden et all 120031; I Chen et all \2004 ) account 
for part of the difference. 

Th e most recent work by Chen et al. ( Chen et all 
120041) describes the process in terms of hard scattering 
from a quark and use Generalized Parton Distributions 
to describe the quark emission and absorption. Chen et 
al. argue that when taking the recoil polarization form 
factors as input, the addition of the two-photon correc- 
tions reproduces i_the. Rosenbluth data. However, Arling- 
ton (jArringtonl . 120051) has shown that when the correc- 
tions of Chen et al. are applied to the new Jefferson Lab 
Rosenbluth data, which have small errors, (see below and 
Figure fT6|) only one-half of the discrepancy is explained. 

Direct tests for the existence of two-photon exchange 
include measurements of the ratio a(e + p)/a(e~p), where 
the real part of the two-photon exchange amplitude leads 
to an enhancement, and in Rosenbluth data where it can 
lead to non-linearities in e. There is no experimental 
evidence of non-lin earities in the Ros enbluth data and the 
e + /e~ ratio data ( Mar et all . Il968t ) are of only modest 



precision, making it difficult to absolutely confirm the 



presence of two-photon effects in these processes. 

It is the imaginary part of the two-photon amplitude 
that can lead to single spin asymmetries but again the 
existing data i Kirkman et all 1 19TCt iPowell et all 1 19T0h 
are of insufficient precision to allow one to make a state- 
ment. There is, however, one observable that has pro- 
vided unambiguous evidence for a two-photo n effect in ep 
elasti c scattering. Groups in both t he US ( Wells et~al\ . 
l200ll ) and Europe (|Maas et aLl . l2005h have measured the 
transverse polarized beam asymmetry. These measure- 
ments are significant but have limited utility in solving 
the G P E discrepancy. The reader interested in more de- 
tail about the existence of two-photon effects and their 
role on the form fact or measurements should refer to 
Ref. (|Arringtonl . [2005l) . 

For quasi-elastic electron-nucleus scattering, the focus 
of this review paper, one therefore must ask the question: 
which form factors should one use in the interpretation 
of the data? What one needs is a parameterization of the 
electron-nucleon cross section for the various kinematics, 
and this parameterization is provided by the form factors 
Ge and Gm determined via Rosenbluth separation of the 
e-N cross sections. Whether the G's from Rosenbluth 
separations or polarization transfer, with or without two- 
photon corrections, are the true (one-photon exchange) 
form factors is largely irrelevant. 

In any case, the discrepancy is largely confined to the 
proton electric form factor G P E . At large momentum 
transfers (above 1 (GcV/c) 2 ) the quasi-elastic cross sec- 
tion is weakly dependent on G E . For example at forward 
angles at 1 (GeV/c) 2 the electric part of the the proton 
cross section is approximately 30%, decreases to less than 
20% at 2 (GeV/c) 2 and to less than 10% at 4 (GeV/c) 2 . 
The electric proton contribution to quasielastic electron- 
nucleus scattering is almost a factor of two smaller due to 
the magnetic contribution of the neutron cross section. 

We also note that another potential problem recently 
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brought up ( Bvstritskiv et all l2007t ). the difference be- 
tween radiative corrections calculated by different au- 
thors, hardly influences the interpretation of quasi-elastic 
electron-nucleus scattering. For the time being, most of 
the data on both the nucleon and nuclei have been ra- 
diatively corrected using the approach of Mo+Tsai (in- 
cluding various improvements made over the years), and 
eventual shortcomings drop out when using nucleon form 
factors in predicting quasi-clastic cross sections. 



VII. SCALING 

In general, the inclusive cross section is a function 
of two independent variables, the momentum transfer 
q and the energy transfer u> of the electron. (The L- 
or T-nature of the scattering might be considered as a 
third independent — discrete — variable) . Scaling refers 
to the dependence of the cross section on a single vari- 
able y(q, u>), itself dependent on q and uj. This scaling 
property is basically a consequence of momentum- and 
energy-conservation in the qu asi- free scattering process 
( Sick et al 1. 119801: IWestl . ll975lh 

Inclusive scattering by a "weakly" interacting probe 
such as the electron can often be interpreted in terms of 
IA (see Sect. HI)) . Quantitative derivations o f scaling in 
IA ha ve been given in several places, see e.g. ( Day et all . 
1990). Here, we first look at a qualitative considera- 
tion, which however contains much of the basic physics. 
Energy and momentum conservation for quasi-free scat- 
tering off an initially bound nucleon with momentum k 
yields 



with 



= [(k + q) 2 + m 2 ] 1/2 - m + E + E, 



ecoil i 



(57) 



with E reco u = lc 2 /(2m(A — 1)). Splitting k into its 
components fcii and k± parallel and perpendicular to 
q and assuming |q|,u> — ► oo, such that the k\ and 
recoil- and removal-energy terms can be neglected, yields 
(uj + m) 2 = fcjj +2&|| |q| + |q| 2 — m 2 . This equation reveals 
that fc|| = y(q, uj), q and u) no longer are independent 
variables. The same value of y (= fen) can result from 
different combinations of |q| and to. The cross section 
cr( |q|,w) divided by the electron-nucleon cross section 
tT eJv(|q|jO-') and a kinematic factor gives a function F(y) 
that only depends on y. This function F has an easy ap- 
proximate interpretation: it represents the probability to 
find in the nucleus a nucleon of momentum component y 
parallel to q. 

The quantitative derivation of ^/-scaling is m ore in- 
volved (|Dav et all Il990l : iPace and Salmd . Il982h . In the 
limit of very large momentum transfer one also finds 
scaling; the main quantitative difference concerns the re- 
striction of the region |k|,i£ which contributes to F(y); 
Fig. \T7\ shows an example. In the limit of very large |q|, 
the scaling function is given by 



(58) 



n(y; |k|) 



dE SQk\,E). 



(59) 



Allowing the upper integration limit to oo, then h would 
correspond to 



i(k) 



dE S(\k\,E) 



(60) 



and F(y) would become the probability n(k\\) to find a 
nucleon with momentum component in the nucleus. 
This approximation is quite reasonable; due to the rapid 
fall-off of S(k, E) with increasing |k| or E, the integral is 
dominated by the region near (|k| = y, E = 0) 



0.4 




F(y) = 2tt / kdk h(y; |k| 



0.4 0.6 
k (GeV/c) 

FIG. 17 Integration region of S(k., E) in Eq. [59] (below corre- 
sponding curves) contribut ing to F(y) at i/=-0.2GeV/c and 
given |q| (|Dav et all Il99(ih . The long-dashed line shows the 
location of the "ridge" of S(k,E) where, at large |k|, most of 
the strength is expected to occur. 



As an illustration of scaling in quasi-elastic electron 
scattering, we show in Fig. [18] some of the available in- 
clusive scattering data for 3 He. The data cover a large 
kinematic range and extend over many orders of mag- 
nitude in cross section, and the quasi-elastic peak shifts 
over a large range of u> with increasing |q|. The same 
data, now shown in terms of the scaling function and 
plotted as a function of the scaling variable y (Fig. [T^]) 
show an impressive scaling behavior for y < 0, i.e. for 
the Iow-lu side of the quasi-elastic peak. The cross sec- 
tion that differs by several orders of magnitude define 
the same function F(y). For y > the values strongly 
diverge, primarily due to the more involved kinematics 
and the different g-dependence of the inelastic e-N cross 
section (A-excitation in particular) which contributes to 
the inclusive cross section at large uj and q. 

In the IA, because of the distribution of strength in E 
of the spectral function, we expect that the asymptotic 
limit of the scaling function would be approached from 
below with increasing |q|. 
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FIG. 18 Cross sections for 3 He(e, e') as function of incident 
electron energy and energy loss, for different incident electron 
energies. 
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FIG. 19 Scaling function for 3 He(e,e'). The various data 
sets are labeled by electron energy (MeV) and angle (deg). 



This scaling property allows the data to be exploited in 
several ways: The presence or absence of scaling tells us 
something about the reaction mechanism (we assumed IA 
without FSI to derive scaling), a residual |q|-dependencc 
of F(y) can tell us something about the |q|-dependence 
of the in-medium nucleon form factor and FSI, and the 
functional form of the experimental scaling function pro- 
vides some insight to the nuclear spectral function. Wc 
address some of these points in more detail below. 

Reaction mechanism. The data show a marked de- 
viation from scaling behavior for y > 0, indicating that 
in this region processes other than quasi-free scattering, 
such as MEC, pion production, A-excitation and DIS 
contribute. This is confirmed by substituting MEC or 
A cross sections for the electron-nucleon cross section in 
our development of F(y). In this instance the data pro- 



foundly fails to scale for y < as well. We thus may 
conclude in particular that, at places where the data ap- 
proximately scale, the contributions of reaction mecha- 
nisms other than quasi-free scattering are smaller than 
the residual non-scaling of the data. Thus, scaling gives 
us direct information on the reaction mechanism, a ques- 
tion we must understand before we can use the data for 
quantitative interpretation. 

Bound-nucleon form factor. The scaling function 
is secured by dividing the experimental cross sections by 
the electron-nucleon cross section. If the |q|-dcpcndencc 
of the electron-nucleon cross section is not correct, the 
data will not scale, at least as long as the range of |q| 
covered is large enough to lead to a large variation of 
the elementary cross section. In order to exploit this 
idea quantitatively, one can compute the scaling func- 
tion using a modified nucleon form factor describing the 
assumed revision of the bound-nucleon form factor, de- 
pending on one parameter, and then fit F(y) using a 
flexible parameterization. The minimum of % 2 of this fit 
gives the best value for the parameter modifying the nu- 
cleon form f actor. It has been found for Iron, for example 
( Sickl . fl986T ) . that within the systematic error of 3% no 
change of the nucleon radius (more precisely the size pa- 
rameter in the Dipole formula) can be supported by the 
data. As the cross sections receive a ~70% contribution 
from magnetic scattering, this result mainly applies to 
the magnetic radius; the limit on any medium modifi- 
cation of the charge radius is twice as large. Given the 
number of models that predict a sizable influence of the 
nuclear medium on the nucleon form factors, the informa- 
tion provided by scaling behavior is quite constraining. 

Constituent mass. In the calculation of the scal- 
ing variable y, the mass m c of the constituent the elec- 
tron scatters from plays an important role, at least as 
long as the recoiling constituent is not ultra-rclativistic 
(recoil energy >3 GeV). The dependence of the scaling 
function on m c can be exploited to learn about the na- 
ture of the constituent. This is important for small uj 
and large q, where it has been suggested that scatter- 
i ng from quark-clusters or individual quarks plays a role 
(|Kumano and Monizl . [l98l iPirner and Varvl . Il98ll ). For 
the kinematical region explored, the data are best ex- 
plained by scattering from nucleons. 

It is only possible to identify the participating con- 
stituent (its mass and form factor) if the data covers 
a large range of momentum transfer. If only a limited 
range is considered, accidental compensations can oc- 
cur which then obscure an interpretation. For instance, 
the simultaneous observation of y-scaling (nucleon mass 
and form factor) and ^-scaling (pointlike, massless con- 
stituents) was rather confusing until it could be shown to 
result from an acc i denta l cancellation of q-de pendencics 
(|Benhar and Liutil . Il995l ; iDav and Sickl . l2004fh 

Nucleon FSI. The final state interaction of the 
knocked out nucleon in general is of minor importance 
in inclusive scattering; the electron carries information 
only about the FSI that takes place within a distance of 
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order l/|q| from the scattering vertex. Subsequent inter- 
actions of the recoil nucleon on its way out of the nucleus 
(which are much more important for e.g. (e, e'p)) do not 
influence the scattered electron. 

At low energy loss FSI does, however, play a role, as 
discussed in Sect. IIIII While the distribution of the spec- 
tral function 5(|k|,i?) in E leads to a convergence of 
F(y, |q|) from below with increasing |q|, the FSI leads 
to a convergence from above. Fig. 1201 which presents 
the scaling function F(y, Q 2 ) for Iron for fixed values 
of y, shows this convergence from above and that, with 
the large momen tum transfers recently made available 
at Jefferson Lab ( Arrington et all . Il999l ). convergence of 
F(y,Q 2 ) can be demonstrated for values of — y as large 
as 0.5 GeV/c. 
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FIG. 20 Convergenc e of F(y,Q 2 ) wi th inc reasing Q 2 , for 
different values of y (|Arrington et all . Il999t ). The rise of 
F(y,Q 2 ) for y = and increasing Q 2 is due to the A- 
contribution. 



kinematics). If at the same time the folding function that 
accounts for FSI becomes q-independent — which is the 
case for NN scattering where the total cross section is 
essentially independent of momentum then the folded 
function also is independent of |q|, and will scale. The 
same observation h as been made for deep inela s tic lep ton- 
nucleon scattering ( Paris and P andharipandc. 20 Q2|). 

It also has been pointed out (jWeinstein and Negeld . 
119821) that for a hard-core interaction the scaling func- 
tion is not directly related to the momentum distribu- 
tion, due to the effect of FSI. However, for less singu- 
lar inte ractions, such as the Paris potential, it has been 
shown ( Butler and Koonfnl Il988| ) via a calculation in 
Brueckner-Goldstone theory that the full response does 
converge to the I A result at large |q|. 

Spectral function at large momentum. The prop- 
erties of S(k, E) at large |k| are closely connected to the 
behavior of the inclusive cross section at large |q| (sev- 
eral GeV/c) and comparatively low uo (several hundred 
McV) . This is qualitatively obvious when considering the 
limit of uj=0. It is only possible (working in the PWIA) 
to transfer a large momentum |q| to the nucleon with 
the result that the nucleon in the final state have both 
small energy, hence small momentum k + q, if the initial 
nucleon had momentum k ~ — q before the scattering. 

The region of the low-w tail of the quasi-elastic re- 
sponse is best studied by considering the ratio of the nu- 
clear and deuteron response (the latter being well known 
experimentally and accurately calculable for any NN po- 
tential). Figure [21] shows the ratio for one of the kinemat- 
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FIG. 21 Per nucleon cross section ratio of nuclear matter and 
the deuteron taken at 3.6GeV and 30°, as a function of the 
scaling variabl e y. The result of the CBF calculation is shown 
as a solid line (|Benhar et all ll994T l . 



Drawing conclusions on FSI starting from scaling 
should be done with caution. Scaling is usually derived 
assuming the IA, so the experimental observation of scal- 
ing seemingly would suggest that FSI are unim portant . 
This, however, is not correct. It has been shown ( Benharl . 
Il999l ) how scaling in the presence of FSI can come about: 
for quasi-elastic scattering the width of the quasi-elastic 
peak becomes constant at large |q| (due to relativistic 



ics at large |q|, small ui where data are available. Data 
and theory are plotted as a function of the variable yn, 
which is basically the component of the nucleon momen- 
tum k parallel to q, calculated for the deuteron. The 
dip at y = results from the fact that the per-nuclcon 
momentum space density at low momenta is higher for 
the deuteron (see Fig. [3]) , the " plateau" in the region of 
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y < —300 McV/c results from the fact that the deuteron 
and nuclear momentum distributions have a similar fall- 
off at large momenta and essentially differ by an overall 
factor. 

Figure [21] shows that the particular nuclear matter 
spect ral function used in this calculation ( Benhar et all 
1 19891) agrees well with the data. Tests carried out by 
renormalizing the spectral function at |k| > kp have 
shown that the cross section ratio at y < —300 MeV/c 
is essentially proportional to >S(|k|,i?) at |k| > kp. One 
needs to realize, however, that for a quantitative study 
of n(|k|) it is important to include FSI, as has been done 
in Fi g. |2~T1 although it was suggested ( Frankfurt et all 
I1993D that the FSI could cancel in this ratio (see also 
Fig. ITTjl . Also, when dealing with nuclei with A>2 it 
is important to use a spectral function S(k, E) rather 
than just a momentum distribution (which ignores the 
independence of S), see Fig. [TT] 

In terms of the scaling variable x the inclusive cross 
sections now reach up to x ~ 3. It is of course tempting 
to interpret the strength near x = 2 (x = 3) as origi- 
nating from scattering off 2 (3)-correlated nucleon sys- 
tems having; 2m (3m ) mass. This type of interpretation 
(|Egivan e£ ^.1 . 120031 ) ■ however, ignores the fact that the 
data exhibit clear j/-scaling, hereby proving that the elec- 
tron scatters from constituents with nucleonic mass and 
nucleonic form factor. The interpretation of cross section 
ratios between nuclei as ratios of correlated strength is 
also hindered by the fact that x, unlike y, is not simply 
related to th e momentum carried by the struck nucleon 
(jLiutil . [l993l ). In addition, the strength at very negative 
y (large x) is strongly affected by A-depcndent FSI, see 
Fig.fJU 

Superscaling. Rec ently, the scaling idea ha s been 
pushed one step further ( Donnelly and Sickl . Il999[ ). Mo- 
tivated by the Fermi-gas model, in which all momentum 
distributions only differ by an overall scale factor — the 
Fermi momentum — Donnelly and Sick have investigated 
whether the scaling functions of different nuclei also can 
be related to each other, by adjusting one overall scale 
factor. It turns out that this is possible for nuclei with 
A > 12; for lighter nuclei deviations near the top of the 
quasi-elastic peak are visible. Fig. [2U shows an example 
the scaling function f(ip') plotted as a function of ip', 
which corresponds to the variable y scaled by a "Fermi 
momentum" . 

This scaling as a function of nuclear mass number ac- 
tually is better realized than the normal scaling which 
is broken by the non-quasi-elastic contributions to the 
response; these non-quasi-elastic contributions are not 
too dissimilar for different nuclei at the same kinematics. 
As the momentum used to scale y is a slow and smooth 
function of A, this superscaling property is particularly 
useful to interpolate data on F(y) in order to predict 
F (y) for nuclei not experimentally investigated. 

For the longitudinal response, superscaling is particu- 
larly well realized, even at large energy loss where the 
transverse quasi-elastic response is obscured by inelastic 
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FIG. 22 Scaling function for nuclei A = 4 - 197 and fixed kine- 
matics (|q| ~lGeV/c) as function oft/)', which corresponds to 
the variable y scaled by a " Fermi momentum" . 



e-N scattering. Fig. [23] shows the longitudinal response 
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FIG. 23 Longitudinal superscaling function extracted 
(|Donnellv and Sickl [i~999n from the available data on C, Ca 
and Fe for momentum transfers between 300 and 570MeV/c. 



function determined via superscaling from the available 
separated data on C, Ca, Fe. This function, incidentally, 
also displays a tail towards larger oj discussed in Sect. ILTl 

This superscaling feature is particularly useful when 
one realizes that (e, e') and ^-induced quasi-elastic pro- 
cesses such as (y, e) or (y, fj,) differ only by the elemen- 
tary vertex, the underlying nuclear physics being the 
same. Cross sections for (e, e') together with superscal- 
ing then allow one to accuratel y predict the cross s ection 
for neutrino- induced reactions ( Amaro et ali . [20051 ). cur- 
rently an area of intense experimental activity. 
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VIII. LIGHT NUCLEI 

Quasi-elastic scattering from light nuclei (A < 4) occu- 
pies a special place. For these systems, several methods 
to calculate the inclusive response have been used, of- 
ten treating both the initial and the final state in less 
approximate ways than imposed by the complexities of 
the many-body system for heavier nuclei. For A=3, for 
example, the first calculation of a realistic spectral func - 
tion was performed 30 years ago ( Dieperink et aLl . ll976t ). 
A variety of approaches to describe quasi-elastic scatter- 
ing from light nuclei is available; we address some of 
their results below. 

The long itudinal response fo r the A=3 nuclei has been 
calculated (|Efros et all l2004f) using the Lorentz integral 



transform technique (see also Sec. HXj) . They use ground 
state wave functions expanded in terms of correlated 
sums of hyper-spherical harmonics, calculated using dif- 
ferent NN forces (e.g. Argonne 18, Bonn- A) and three- 
body force (3BF) models. As their calculation is essen- 
tially nonrclativistic, they restrict it to |q| <500 MeV/c, 
and consider the longitudinal response to avoid the diffi- 
culty of MEC. Fig. |2"41 shows their results for 3 H and 3 He 
compared to data. 
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FIG. 24 Compariso n of experimental and theoretical response 
(|Efros et all I2004T ) for Argonne 18 potential (dotted) and 
Argonne V18+3B F (solid). Da ta are shown as open circles 
(|Dow et afill988f ) and squares (|Marchand et af.l Il985h . 



This calculation finds rather small differences between 
the responses calculated using different modern NN 
forces. The three-body force leads to a systematic re- 
duction in the height of the quasi-elastic peak, presum- 
ably as a consequence of the tighter binding which leads 
to a more compressed (extended) distribution in radial 
(momcntum)-space. While for 3 He the 3BF is helpful 
in explaining the data, this does not seem t o be the case 
for th e 3 H nucleus. In a very recent paper ( Efros et all 
I2005D the range of applicability of this nonrclativistic cal- 
culation has been extended. 

The inclusive cross section has also been calcu- 
lated separa t ely fo r the 2-body and 3-body breakup 
dGolak et ad 1 19951) . 7or the ground state, Golak et al. 



use a solution of the 34-channcl Faddecv equation for 
the Bonn-B NN-potential. For the final state, the au- 
thors separate the contribution from the (symmetrized) 
plane-wave approach and the one from rescattering pro- 
cesses, summed to all orders. This calculation again is 
non-relativistic, and does not include MEC. 

Fig. [231 shows the results of Golak et al. for 3 He and 
3 H , with the contributions of 2- and 3-body breakup sep- 
arately displayed. 
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FIG. 25 Comparison of theoretical (|Golak et all 1 19951 ) and 
experimental response at — q — =300MeV/c for 3 i/e (a) and 
Z H (b). Dotted: 2-body breakup, da shed: thre e -bod 
breakup. Da ta are shown as open circles (|Dow et all 1 198 ~ 
and squares (|Marchand et al.lll985h . 



The inclusive cross section has also been calculated 
using a realistic coupled-channel potential with single A- 
isobar excitation for the initial and final hadronic states, 
with the corr esponding e.m . curr ent with two-baryon 
contributions i Deltuva et aLl . [2004h . The potential is an 
extension of the CD-Bonn NN-potential. The main A- 
isobar effects on observables then result from the effective 
three-nucleon force and the corresponding effective two- 
and three-nucleon exchange currents. Both initial bound 
state and final continuum state are exact solutions of the 
3-particle scattering equations. 

In particular at low momentum transfer Deltuva et 
al. find large effects of the A in the threshold region of 
the transverse response. At the larger momentum trans- 
fers, e.g. at 500 MeV/c shown in Fig.[2S]for 3 He, the con- 
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tribution of A degrees of freedom is smaller, in agreement 
with the observation made in Sec. lIXl The small shift be- 
tween data and calculation has been assigned to the use 
of non-relativistic kinematics (necessary for consistency) . 
Deltuva et al. also have calculated other e.m. observ- 
ables such as elastic form factors and exclusive quasi- 
clastic cross sections, and find rather large changes when 
allowing for the coupled channel nuclcon-A case. 
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FIG. 26 Comparison of theoretical (jDeltuva et al 
and experimental response at |q|=500MeV/c. Solid: 
coupled-channel nucleon-A, dashed: purely n ucleonic case. 
Data are shown as ci rcles (|Dow et al. |, 1 1988ft and squares 
(jMarchand et all Il985l ) . 



As pointed out in Sec. IVIII the scaling function F(y) is 
closely related to the spectral function and momentum 
distribution. In the limit of very large momentum trans- 
fer, and when neglecting the effects due to finite E and 
FSI, the scaling function is given by 



(61) 



F(y) = 2n n(|k|)d|k| 
J\y\ 



in which case the momentum distribution can be ob- 
tained from 

1 dF 
2iry dy ' 

For the deuteron, neglecting the effects of finite E and 
FSI is not unreasonable as the data reach very high q. 



n(\k\) 



with k 



y 



(62) 



The n(k) derived via scaling (jCiofi degli Atti et alj . ll987t ) 
is compared in Fig. [57|to the ones obtained from (e, e'p) 
reactions and theory. 

For heavier nuclei, both FSI and the effects of the dis- 
tribution in E are no longer negligible, and can only 
be incorporated in predictions of F(y, Q 2 ) starting from 
S(k,E). 




FIG. 27 Comparison of deuteron momentum distribution ob- 
tained from (e,e') (squares), ( e,e'p) (triangles) and theory 



(curves) i|Ciofi degli Atti et all \W8 



IX. EUCLIDEAN RESPONSE 

In previous sections, we have discussed quasi-elastic 
scattering in terms of PWIA, with corrections for FSI. 
This type of theoretical description often is the only prac- 
tical one, as a less approximate treatment of the final, 
relativistic, continuum state is hard to come by. 

For comparatively low momentum transfers, alterna- 
tive and more reliable treatments are possible. For A = 3 
and A = oo th e response in the non-relativistic region can 
be calculated ( Fabrocini and Fantonil . Il989t iGolak et all 
I1995D . For nuclei in between, integrals over the quasi- 
clastic response can be studied. These integrals can be 
expressed as expectation values of the ground-state wave 
function, without need for an explicit treatment of the 
continuum state. This approach has been followed via 
both the Lorentz transf orm technique ( Efros et q/J . |l994l; 
iLeidemann et al. 1 19971). and t he us e of the Laplace trans- 
form ( Carlson and Schiavillal . [l992l) . The most extensive 
results, including one- and t wo-body currents, ar e avail- 
able for the latter approach ( Carlson et all |2002[ ) which 
we discuss in more detail below. 

The Euclidean response is defined as an integral over 
the quasi-elastic response 



E 



T,L 



|q|,r) 



exp[-wr] R TL Qq\,u) da; (63) 



with t = Q 2 /Am. In the Lorentz transform technique 
an additional Lorcntzian factor which enhances the in- 
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tcgrand at a given u> is used. The Lorentz transform 
technique has the advantage that the response as a func- 
tion of lj can be reconstructed by the inverse transfor- 
mation, something that is not practical for the Euclidean 
response. 

The longitudinal and transverse Euclidean response 
functions represent weighted sums of the corresponding 
i?i(|q|,tj) and Rt{\(\\ 1 uj): at t=0 they correspond to the 
Coulomb and transverse sum rules, respectively, while 
their derivatives with respect to r evaluated at t=0 cor- 
respond to the energy- weighted sum rules. Larger values 
of r correspond to integrals over progressively lower en- 
ergy loss regions of the response. 

In a non-relativistic picture, the Et,l can be simply 
obtained from: 

E L (\q\,r) = (0|pt(q) exp [-(F-.Eh)r]p(q)|0) 



0.020 



exp 



q 2 r 
2 Am 



|(0(q)|p(q)|0)| 2 , (64) 



and similarly for £r(|q|, t), with the charge operator 
p(q) replaced by the current operator jr(q)- The elastic 
contributions have been explicitly subtracted, and |0(q)) 
represents the ground state recoiling with momentum q. 

The study of the Euclidean response has the outstand- 
ing advantage that _E(|q|,r) can be calculated from the 
ground state properties alone; no explicit treatment of 
the final continuum state is required. For the A=3,4 
ground states, very precise wave functions are available, 
and the effects of MEC can be included using the two- 
body operators well established in elastic and inelastic 
electron scattering from light nuclei (for a review sec Sick, 
2001 ). 

The Euclidean response has the disadvantage that we 
lack an intuitive in terpretation of this i ntegrated quan- 
tity. Model studies (jCarlson et all l2002f) show the sensi- 
tivity of E(\q\ 7 r) to properties of i?(|q|, w), see Fig. [28] 
The top figure shows various modifications of the 'nor- 
mal' response, the bottom figure shows the effect upon 
the ratios of the resulting Euclidean responses to the 'nor- 
mal' one. These studies show that, for the responses that 
can be extracted from the data, the region 0.01 < r < 
0.05 MeV -1 is the most relevant one for a comparison 
with theory. Below r = 0.01 MeV -1 the contribution 
of the tail of the A in the experimental response is too 
important, above r = 0.05 MeV -1 the response is totally 
dominated by the contribution of v ery small values of lo . 

The calculations discussed here ( Carlson et all l2002h 
have used the standard expressions for the one-body 
electromagnetic operators, obtained from a rclativis- 
tic reduction of the covariant single- nuclcon current. 
The two-body current operator consists of "model- 
independent" and "model-depe ndent" comp onents, in the 
standard classification scheme dRiskal.ll989h. The model- 
independent terms are obtained ( Schiavilla et all 1 1990h 
from the nucleon-nucleon interaction. For the model- 
dependent pieces, the calculation includes the isoscalar 
p7T7 and isovector ujttj transition currents as well as the 
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FIG. 28 Model responses derived from the "normal" response 
plus various modifications (top) . The bottom figure gives the 
corresponding Euclidean response in terms of the ratio to the 
one from the "normal" one. 



isovector current associated with excitation of interme- 
diate A-isobar resonances. The two-body charge oper- 
ators ( Schiavilla et all Il990h include the 7T-, p- : and uj- 
mcson exchange charge operators, the (isoscalar) pix^j and 
(isovector) uj-k^i couplings and the single-nucleon Darwin- 
Foldy and spin-orbit relativistic corrections. 

We show in Fig. [2U the results obtained for 3 He and 
He at one value of |q|. The ground-state wave functions 
used in this stu dy were obtained wit h Variational Monte 
Carlo approach ( Carlson et all [20021 ) and the Argonne v$ 
N-N interaction plus the UIX three-nucleon interaction. 

The helium nuclei studied by Carlson et al. are of 
particular interest as for A He the excess in the transverse 
strength is maximal among all nuclei, and grows by a 
factor of two between A=3 and A=4. This excess - 
presumably due to MEC — had not been understood in 
the past; the many calculations of MEC for a multitude 
of nuclei gave results that were rather discordant, and 
always much too small. The data for the Helium nuclei 
also show, that this excess covers the entire quasi-elastic 
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peak, and not only the "dip" region between the quasi- 
elastic- and A-peak. 
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FIG. 29 Transverse Euclidean response divided by the proton 
magnetic form factor for 3 He (top) and 4 He at 600 MeV/c 
momentum transfer, scaled by e TUJq " . Data (+), IA (o) and 
full calculation (x). 



Fig- HS1 shows that the full calculation, which includes 
MEC, is in good agreement with the data (we pointed 
out above why the region r < 0.01 McV -1 should be 
ignored). The calculation predicts rather accurately the 
enhancement of the transverse strength due to MEC and 
the doubling between A=3,4, and it also does quite well 
in predicting the <5 2 -dependence (not shown). 

The good reproduction of the transverse strength 
at first sight comes somewhat as a surprise, 
given the lack of suc c ess o f previous MEC- 
calcu lat ions (lAlberico et all 1984 lAmaro and Lallena l. 



1992 



1989 



Angui ano et all 119961 ; iBlunden and ButlerL 



Carlson and Schiavilla 



Dckker et all 



1994 


Donnelly et al. 


19781; iFabrocini, 


1997^ 


Gadivak and Dmitriev. 


1998 


; Kohno and Ohtsuka 


, 1981 


Leiclcmann and Orlandini, 


1990i; Van der Sluvs et al\ 



19951 ; IVan Orden and Donnelly! , Il98ll ). Carlson et 



al. have therefore investigated more in detail the 
reason for the large MEC contribution. They have 
found that, in agreement with pre vious stud i es that 
have received too little attent ion ( Fabrociml Il997bt 
iLeidemann and Orlandinil . Il990t ). MEC only produce 
large effects in combination with ground-state wave 
functions calculated including the short-range n-p cor- 
relations. As most previous calculations were based on 
independent-particle type wave functions, the smallncss 
of the resulting MEC contributions thus is understood. 
To verify this point further, Carlson et al. have repeated 
their calculation using the same operators, but with 
a Fermi-gas wave function. Instead of a enhancement 
factor of 1.47 coming from MEC at |q| =600 McV/c 
they find a factor of 1.06 only, i.e. an eight times smaller 
MEC effect. 

The results of Carlson et al. also show, somewhat sur- 
prisingly, that the MEC contribution is large at low mo- 
mentum transfer. It decreases towards the larger Q 2 , in 
agreem ent with the ex pectation that at very large Q 2 
it falls ( Sargsiarl l200ll ) like Q~ 4 relative to quasi-elastic 
scattering. 

From the above discussion it becomes clear that the 
Euclidean response, despite inherent drawbacks, is a very 
valuable quantity. Since the final continuum state does 
not have to be treated explicitly, calculations of much 
higher quality can be performed than for the response, 
and the role of two-body currents can be treated quan- 
titatively. The comparison between data and calculation 
has shown in particular that for a successful prediction 
of MEC correlated wave functions for the ground state 
are needed; such wave functions today are available up 
to A ~ 12 and for A = oo. Unfortunately, the usage of 
the Euclidean response for the time being is restricted to 
a regime where relativistic effects arc not too large, such 
that they can be included as corrections. 



X. L/T-SEPARATION AND COULOMB SUM-RULE 

In impulse approximation, and when neglecting the 
(small) contribution from nucleonic convection currents, 
the longitudinal and transverse response functions Rl 
and Rt contain the same information and have the same 
size. This has sometimes been called "scaling of the 
th kind" (see Sec lVII[) . It has been realized early on, 
however, that the transverse response receives signifi- 
cant contributions from meson exchange currents and 
A-excitation (which are of largely transverse nature). It 
therefore is clear that there is a high premium on separat- 
ing the L- and T-responses, both because the L-response 
is easier to interpret and because of the additional infor- 
mation contained in the T-response. 

The separation of the L- and T-responses is performed 
using the Rosenbluth technique, which is justified only 
in the single-photon exchange approximation. The cross 
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FIG. 30 Longit udinal (low e r dat a set) and transverse re- 
sponses of 12 C (|Finn et al\ , I1984T ), plotted in terms of the 
scaling function F(y). 



section, divided by a number of kincmatical factors 

J^r ^— % = e WW,") + MM,*) = s 

diirfw <7 M ott Q 



(65) 



is a linear function of the virtual photon polarization 



e=(l + -^-taa»-) 



(66) 



with q (Q) being the 3-(4-)momentum transfer and e 
varying from to 1 for scattering angles 9 between 180° 
and 0°. The slope of the linear function yields R^, the 
intercept at e = yields Rt- Fig. [3D] shows an early 
example for an L/T-separation, and demonstrates the 
excess observed for the transverse strength. 

While conceptually very straightforward, this L/T- 
separation is difficult in practice. It involves data tak- 
ing at the same |q|, but varying e, i.e. varying beam 
energy. For an accurate separation of Rl and Rt obvi- 
ously the largest possible range in e, hence beam energy, 
is required. As data are usually not taken at constant 
|q|, but at a given beam energy and variable energy loss, 
obtaining the responses at constant |q| involves interpo- 
lations of the data. We show in Fig. [31] two examples 
for a Roscnbluth separation, performed on the low- and 
large-w side of the quasi-elastic peak, which also illus- 
trate the importance of the forward angle (high energy) 
data for the determination of Rl, i.e. the slope of the 
fit. 

The Rosenbluth technique is applicable in Plane Wave 
Born Approximation PWBA, and fails once Coulomb dis- 
tortion of the electron waves is present. Neglect of distor- 
tion is justified for the lightest nuclei alone, and only if 
Rt is not much bigger (or much smaller) than Rl . When 
one of the two contributions gets too small, even mi- 
nor corrections due to Coulomb distortion can have large 
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FIG. 31 Rosenb luth separation for 56 Fe and |q| =570 MeV/c 
(|JourdarJ . ll996l ). The dashed lines are fits to the Saclay data 
alone, the solid lines are fits to the world data which include 
the forward-angle SLAC data and emphasizes the importance 
of a large range in e. 



effects. At large |q|, for instance, even the determina- 
tion of the proton charge form factor via the Rosenbluth 
technique is significantly aff ected by Coulomb corrections 
( Arrington and Sickl . 12004 ) . In order to extract Rl and 
Rt in the presence of Coulomb distortion, the data have 
first to be corrected for these effects; this is discussed in 
Sec. KB 

Here, we concentrate on the discussion of the longi- 
tudinal response. For practical reasons, the determina- 
tion of the longitudinal response is possible only in a 
g-range that is somewhat limited. In addition, at low 
|q|, typically below twice the Fermi momentum, the re- 
sponse is affected by Pauli blocking, which in many of 
the approaches used to describe inclusive scattering is 
not properly treated. At large |q|, typically above 0.8 
GeV/c, the transverse response dominates the cross sec- 



tion, both due to the |q| 2 /Q 2 factor in Eq. (f6"5"l) and 
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the increasing A contribution (see Fig. [30]) . such that 
an accurate determination of Rl becomes very difficult. 
Energy-dependent experimental systematic errors must 
be handled with great care. 

One particular use of the longitudinal response has re- 
ceived much attention: the determination of the Coulomb 
Sum Rule (CSR). In the non-relativistic regime, when 
short-range correlations between nucleons and the effect 
of Pauli blocking is neglected, the CSR takes the simple 
form 



Sl{H) 



G e 2 



du> 



(67) 



where G e 



is the threshold for 



(G%+G 2 en N/Z) and ^ 
particle emission. In the limit of large |q| Sl should be 
one. To say it in words: when neglecting the small con- 
tribution from the neutron charge form factor G en , the 
integral over the longitudinal response counts the num- 
ber of protons times the square of the proton charge form 
factor G\ y . 

The history of the CSR is a very checkered one. Early 
work proposed the CSR as a tool to study short-range 
correlations SRC betw een nucleons ( Czvz and Gottfried! . 



119631 : ICottfriedl . ll963D . These correlations move strength 
to large energy loss, and partly out of the physical re- 
gion; at very large |q|, these correlation contributions go 
to zero in the sum r ule. The series of L/T -separations 
performed at Bates ( Altemus et all Il980f) and Saclay 



( Meziani et a/.l . ll984j ) found effects in Sl that were much 



bigger than could be expected from SRC: in the region 
of |q|=350-550 MeV/c S L was up to 50% lower than 
expected, the deficit increasing with increasing |q| and 
increasing nuclear mass number A. These observations 
have widely been interpreted as a medium-modification 
of the proton charge form factor ( Celenza et all \l98& 
iMuldersl . Il986l : iNobld . Il98lh . 

The experiments dealig with ther CSR have absorbed 
much of the attenti on. It was not g enera lly known that 
some experiments ( Altemus et all Il980i ) had suffered 



from rescattcring of electrons on the "sn out" connect 



ing s cattering chamber and spectrometer ijDeadv et 



19861 ) (see also Sec . IIVI). It also h as taken a long time 



before a reanalysis ( Jourdanl . Il996l ) removed a number of 



deficiencies in the Saclay analysis of the data (see below). 

L/T separations have been have been performed 
for a number of nuclei and momentum trans- 



fers. (lAltemus et all 19801 Barreau et all 


1983 


Blatchlev et all 19861 Chen et all 1991; Deadv et ai. 


19861. 19831: Finn et ai, 1984: iMeziani et all 


1992 


VanReden et ai. 


19901: IWilliamson et all 


1997 


Yates et ai, 19931: 


Zimmerman!. Il969l). Often 


(but 



not always, see e.g. Williamson, 1997), these separations 
were quite limited in the e-range as data from one 
facility only were included; from the determination of 
the proton form factors G ep and G mp it is well known 
that the more reliable results are obtained from an 
analysis of the world data spanning the largest possible 
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FIG. 32 Separated response functions for 56 Fe and 
|q|=570 MeV/c. The solid curve corres ponds to the CBF cal- 
culation ([Fabrocini and Fantonil . 1 19891 ) . 



range in e (electron energy). In almost all cases, only 
approximate Coulomb corrections were included, using 
the EMA version of the effective-momentum transfer 
approach (see Sec. IXI[) . Particularly for data producing 
low Sl the longitudinal response functions show an 
unphysical behavior at large uj: they dive steeply 
towards zero, and would, for responses without discon- 
tinuities, obviously be negative just beyond the range 
of uj shown (s ee e.g. Jourdan ,1 996). From microscopic 
calcu l ations (|Dellafiore et all 19851 : iDieperink et all 
119761 : iFabrocini and Fantonil . Il989h using realistic nu- 
clear spectral functions we do know, however, that the 
response at large uj should approach zero slowly, much 
more slowly than the response at low uj, as a consequence 
of the components of large removal energy E present in 
realistic spectral functions S(k, E) (see e.g. Fig. [24)) . 



The analysis of th e data performed by (jCarlson et all 
120031 : iJourdanl . Il996f ) included the world data to employ 
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the largest e-range and improved a number of aspects of 
the Saclay analysis. In particular, the integral in Eq. (|67p 
was divided by the correct G ep , the well-known rela- 



tivistic corrections were included (de Forest! . fl98l . and 
Coulomb corrections in DWBA were made; in addition, 
the contribution of Rl(lu) above the upper integration 
limit — which experimentally is far from oo due to the 
limited range of data on S'l(w) — was added. These 
L/T-separations extend to |q| ~ 600 MeV/c, only for 
A He, where the overlap of quasi-elastic and A strength is 
a lesser problem, the L-strength is known up to 1 GeV/c. 
The corrections to the 1984 analysis, which all happen to 
go in the same direction, increase the CSR by a substan- 
tial amount. For example, the CSR for the largest |q| 
and heav i est nu cleus, 56 Fe and 570 MeV/c analyzed in 
(jjourdanl . [l99ll . amounts to 0.98±0.15. Fig. [521 shows 
the corresponding L- and T-responses. 

In a subsequent reanalysis of the (e,e') data 
( Morgenstern and Mezianl . 1200 ll ) resulted in a less ex- 
treme reduction in the CSR result, for Iron the CSR was 
found to be 82% of the expected value. This result for 
the sumrule is still significantly smaller than the one of 
( Jourdar] . Il996l ). the difference being largely due to the 
use of EMA for the Coulomb corrections. This empha- 
sizes that the Coulomb corrections are very important 
for the determination of the CSR; particularly for the 
largest q-valu e of 570MeV/c and the large uj, where the 
differ ence of ( Jourdanl . Il997t iMorgenstern and Mezianil . 
2001) is largest, the backward angle data go down to 
scattered-electron energies a s low as 130MeV, where ap- 
proximations like EMA fail (jAste et aLl . 120051 ). 



The issue of medium modifications of proton electro- 
magnetic structure has been recently revived by the re- 
sults of a polarization transfer 4 He(e, e'p) 3 H measure - 
ment carried out at Jefferson Lab ( Strauch et al I l2003h . 
However, the interpretation of the experimental data in 
terms of medium-modified form factors is challenged by 
the results of a theoretical calculation, carried out us- 
ing accurate three- and four-nucleon bound-state wave 
functions, a realistic model for the nuclear electromag- 
netic current operator and a treatme nt of final-state- 
inter actions with an optical potential ((Schiavilla et all 
I2005D . In Rcf. (|Schiavilla et all |2005j no significant 
discrepancies are found between theory and experiment, 
both for the ratio of transverse to longitudinal polar- 
ization transfers and for the induced polarization, when 
free-nucleon electromagnetic form factors are used in the 
current operator. 

To close the discussion of the Coulomb sum, we add 
as a caveat that even for the longitudinal response the 
contribution of MEC is not entirely neg ligible. In the 
_[c(J_~ 1 GeV/c region it has been shown (jCarlson et 



120031 ) that for 4 He the contributions are of the order of 
10%. 



XI. COULOMB CORRECTIONS 

The effects of the static Coulomb field of nuclei upon 
quasi-elastic scattering has posed a persistent problem. 
The presence of the Coulomb potential — for Lead of 
order of 25-30 MeV in the nuclear interior — has a ma- 
jor effect on quasi-elastic cross section measured in the 
several-hundred MeV energy region. It invalidates the 
linear relation (|65|) used to separate the longitudinal and 
transverse responses using the Rosenbluth technique. 

While the treatment of the Coulomb distortion of the 
electron waves in the quasi-elastic region presents no con- 
ceptual problems, the practical application has been dif- 
ficult. Reliable calculations of Coulomb distortion have 
not been easily accessible to analyses of experimental 
data. As a consequence, most experiments have been 
analyzed without considering Coulomb distortion effects, 
or by using relatively simple recipes. 

The Coulomb distortion is more important for the 
lower-energy data and the heavier nuclei. Its effects are 
visible in the separations of longitudinal and transverse 
strength, where the small contribution — in general the 
longitudinal one — is most affected. It is likely that one 
of the main problems with the longitudinal strength — 
the "diving to negative res ponses" mentioned in Sec. [X] 
— is related to this aspect ( Traini et all 1 1988h . 

The Coulomb distortion can be treated in Distorted 
Wave Born Approximation (DWBA) by using electron 
wave functions calculated as solutions of the Dirac equa- 
tion for the known nuclear charge distribution; differ- 
ent p r ograms ha v e been employed (ICo and Heisenbergj . 
119871 ; IKim et all Il996l lUdias et all Il993l). m_partic- 
ular by the group of Onley and W right ( Jin et all , 
Il992l ; IZamani-Noor and Onlevl . Il986l ). In the s e cal- 
culations, high-energy approximations (iKnolll . Il974t 
iLenz and Rosenfelderl Il97lt iRosenfelderj . Il98dl ) can be 
employed, as in general the electron energy is much 
higher than the Coulomb potential. Such calculations 
have been successfully carried out, and can serve as a 
benchmark. For a systematic analysis of experimental 
data, the numerical effort is usually too big and un- 
practical. This is true in particular as often computer 
codes have been developed for (e,e'p), in which case the 
(e, e') cross section has to be generated by summing over 
all possible initial and final states of the knocked-out nu- 
cleon. 

Simplifications are possible when going t o the Eikona l 
approximation dGiusti and Pacatl . Il987l ; iTraini et all 
119881 : lYennie et all Il965h . When using the lowest or- 
der expansion in Za, one finds two dominating effects: 

1. As a consequence of the attractive electron- nucleus 
Coulomb interaction, the effective energy of the incident 
and scattered electron at the moment of the scattering 
is increased by the Coulomb potential V. This has the 
consequence that the effective momentum transfer q e ff 
squared, on which the response functions R depend, is 
increased by a factor (E e + V){E e ' + V)/E e E e >. 

2. Due to the attractive Coulomb interaction, the clcc- 
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tron "plane" waves are focused onto the nucleus, hereby 
increasing the wave function at the location of the nu- 
cleus, with a corresponding increase in the scattering 
cross section. 

In the calculation of the inclusive cross section, the 
overall effect of these corrections is accounted for by us- 
ing the Mott cross section with the unmodified electron 
energy, but with q e ff as an argument of the response 
function. This approximation has been baptized Effec- 
tive Momentum Approximation (EMA). Versions in the 
literature differ by the choice of V. The results of EMA 
have been compared to the ones fr om DWBA calculations 
(iJourdanl . Il997l : iKim et ad . fl996h . 



Hig her order terms have been included ijTraini et all , 
Il988h . These terms have been calculated using severe 
approximations in the expansio n around r = 0, and are 
not recommended ( Trainil . l200l[ ) . Due to the approxima- 
tions, the second-order effects found are about as large 
as the first order, which is indicative of problems of the 
expansion. 

The term " EMA" is often confused in the literature be- 
cause of the fact that it is used for two different choices 
of the nuclear Coulomb potential V. This parameter V 
is often evaluated for the nuclear center r = assuming a 
homogeneous nuclear charge density, in which case V = 
Vq = ZZot/2R eq , with R eq being the equivalent radius 
(a good ap p roxim ation being [1.1A 1 / 3 + 0.86A _1 / 3 ]/m 
( Kim et all , Il996l )'). It has been recognized early on, 
however, that a better choi ce would be an ap propriate 
average Coulomb potential I Rosenfelderl . 1 19801 ). For this 
reason, many applications of EMA use for V the value at 
the nuclear surface, V s = (2/3) Vo, where most nuclcons 
are located. 

In order to improve upon the quality of the Coulomb 
corrections without resorting to the full solution of the 
Dirac equation, one can employ the Eikonal Distorted 
Wave Born Approximation eDWBA, where the electron 
waves, in the DWBA approach solutions of the Dirac 
equa tion, are calculat ed using the Eikonal approxima- 
tion ( Aste et oil . l2004f ) . In this case the electron current 
is modified by an additional Eikonal phase and a change 
in amplitude. This type of calculation can be carried out 
for realistic shapes of the nuclear Coulomb potential, and 
eDWBA, contrary to the full solutions of the Dirac equa- 
tion, can be extended more easily to the larger energies 
of interest for modern experiments. 

While eDWBA is a fairly practical approach that can 
be employed on a routine basis, it still is much more in- 
volved than EMA. Aste and Jourdan have identified one 
problem of EMA, an d have introduced an EM A-like ap- 
proach called EMA' (jAste and Jourdanl . 12003) . In EMA, 
one tries to treat two distinct effects mentioned above: 
the increase of the electron momentum due to Coulomb 
interaction with the nucleus, and the focusing of the elec- 
tron waves. To handle the increase of the electron mo- 
mentum (and momentum transfer) it clearly makes sense 
to use a potential V that corresponds to the Coulomb 
potential averaged over all nucleons; here the use of the 



potential at the nuclear surface is a good approximation. 
For the focusing effect, on the other side, the value for 
the nuclear center is a better approximation, as the fo- 
cusing takes place all along the trajectory of the electron 
through the nucleus. As the electron approaches (leaves) 
the nuclear center the focusing is smaller (larger) . Hence 
the value at the center is a good compromise. 

Subsequent studies using solutions of the Dirac equa- 
tion in the Coulomb field of the nucleus ( Aste etail , 
I2005D indicate that, in addition to the enhanced focus- 
ing in the longitudinal direction accounted for by EMA', 
there is a reduction for non-central electron trajectories. 
The two effects together give an overall focusing which is 
not far from the one obtained with EMA. 

Recent studie s by Tjon and Wallace 

( Tion and Wallace! l2006f ). performed using the Eikonal 
expansion for the electron wave function, indicate that 
EMA somewhat overestimates the Coulomb effects. The 
authors give a recipe for correcting the deficiencies, 
which however is not very practical for analys es of data. 

The recent work of (|Kim and Wright! . l2005h . who per- 
formed calculations of (e,e) in both full DWBA and also 
in an approximate version DW which was found to agree 
well with DWBA, indicates that EMA is a good approx- 
imation for the transverse part, but that it works poorly 
for a longitudinal cross section. The authors, however, 
offer no explanation for the difference. 

It would clearly be desirable to certify, via exact cal- 
culations, the validity of some EM A- type approach, as 
only such an approach could also be applied to the im- 
portant region of the large-w-side of the quasi-elastic 
peak, where one has to deal with the overlap with a large 
A-contribution, the Coulomb corrections for which also 
need to be dealt with properly. 

The ideal way to experimentally check the Coulomb 
corrections is a comparison of electron and positron scat- 
tering. Unfortunately, positron beams are hard to come 
by, and experiments with the secondary positron beams 
are much more difficult than wit h electrons. One su ch ex- 
periment has been carried out ( Gueve et all [l999) . The 
data, unfortunate ly, suffer from normalization problems 
( Aste et al\ , 12004 ) ; from the position of the quasi-elastic 
peak one can, however, deduce that it is appropriate to 
use V s for the calculation of q e f / . 



XII. NUCLEAR MATTER 

For nuclear matter the Schrodinger equation for nu- 
clcons bound by the nucleon - nucleon interaction (de- 
duced from NN scattering) can be solved with very few 
approximations. The translationally invariant nature of 
the medium, where solutions can all be written in terms 
of plane waves, simplifies the calculation very much. As a 
consequence, the quality of nuclear matter wave functions 
is comparable to the one for the A=2,3,4-nuclei. Due to 
the " exact" nature of the nuclear matter single particle 
wave function, both the long-range and short range prop- 
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erties are well under control. This is in contrast to finite 
nuclei, where calculations, that are designed to do well on 
the long-range properties (mean-field calculations), usu- 
ally do badly on the short-range aspects. 

Unfortunately, only some integral properties such as 
density and binding energy are known experimentally; 
little is known on the short-range properties of nuclear 
matter. Quasi-clastic scattering here provides valuable 
information. 

It is important to realize that inclusive electron scat- 
tering at large Q 2 is sensitive only to rather " local" prop- 
erties of the medium. The spatial resolution of (e, e') is of 
order l/|q|, which at large transfer is small. In particular, 
the scattered electron is not sensitive to the interactions 
of the recoiling nucleon outside this range. This allows 
for an extrapolation from finite-nucleus data using the 
LDA. 

The extrapolation procedure ( Day et all . 1 1989T ) starts 
from the consideration that the nuclear response is essen- 
tially the incoherent sum of contributions from individual 
nucleons. As the average value of the density in the nu- 
clear interior and the shape of the density distribution 
in the nuclear surface are approximately /l-indcpendcnt, 
the response can be divided into a volume component, 
proportional to the mass number A, and a nuclear sur- 
face component proportional to A 2 / 3 . It is the former 
one that is of interest when discussing nuclear matter. 
The ratio of the surface to volume contributions is thus 
proportional to A -1 / 3 . Extrapolation of the nuclear re- 
sponse per nucleon to A^ 1 ^ = {A — > oo) as a linear 
function of A -1 / 3 yields the nuclear matter response. 

In order to illustrate this approach we reproduce in 
Fig. [33J o ne example for a n extrapolation as a function 
of A- 1 / 3 (|Dav et al.l . ll989l ). Ignoring 4 He (for which the 
properties of the density cited above are not valid) , the 
nuclear response s(q 7 uj) for A = 12-197 is well fit by 
a linear function of A -1 / 3 . The plot on the bottom of 
Fig. [33] gives the same extrapolation as a function of A. 
This figure reveals that the extrapolation as a function 
of A is unwieldy though the curve better imparts the sat- 
uration of the response. Even heavy nuclei significantly 
differ from nuclear matter due to the large fraction of 
surface-nucleons . 

In order to obtain the response function for symmetric 
nuclear matter, one makes use of additional knowledge. 
For quasielastic scattering, the relative contribution of 
protons and neutrons changes as a function of A. Al- 
though the protons dominate due to the larger electron- 
proton cross section, the contribution of neutrons is not 
negligible. For the extrapolation, one assumes that the 
response functions for protons and neutrons are the same, 
and the trivial dependence on N,Z is removed by extrap- 
olating the quantity 



s(q, uj) = a(q, u) / (Za ep + Na en ) 



(68) 




3 
& 



FIG. 33 Extrapolation of nuclear response per nucleon at 
fixed q and u (E = 3.6 GeV, 6 = 16°, and u = 180 MeV) 
as a function of J 4~ 1/ ' 3 (top) and as a function of A (bottom) 
where the extrapolated value of the nuclear matter response 
is indicated by the arrow. 



This was made possible by the availa bility of data for 
nuclei with A = 4, 12, 27, 56 and 197 (|Dav et all Il993l ) 
taken at the same incident energies and angles. Recent 
data taken at Jefferson Lab with A = 3, 4, 9, 12, 63 and 
197 will allow this to be extended over an even larger 
range of q and to. 



The nuclear ma tter response has been extrapolated 
( Day et a71 . ll989l ) from the nuclear response measured for 
finite nuclei over a large region of q and w, see Fig. 1341 



XIII. RELATED AREAS 

Inclusive scattering from co mposite systems i s use d 
as a tool in a number of areas (|Silver and Sokoll . Il989h . 
The corresponding processes have many aspects that are 
closely related to quasi-elastic scattering from nuclei, but 
they also exhibit significant differences. We now address 
those areas where inclusive scattering has been harnessed 
to study diverse composite systems. 

These areas differ not only by the nature of the com- 
posite system investigated, but also by the probe used: 
photons, low-energy electrons, neutrons, high energy 
electrons, muons and neutrinos. The energies of the 
probes cover many orders of magnitude, from mcV to 
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FIG. 34 Nuclea r matter response at momentum transfers up 
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mentum of 0.5 a.u. is due to the free-electron Fermi gas 
and the tail at larger momenta is du e to inter-electron- 
interactions ( Eisenberger et adll972l ). 
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GeV, scaling with the dimensions relevant for the com- 
posite systems. 

Historically, the first area where inclusive scattering 
became prominent was the measurement of the Compton 
profile, i.e. quasi-elastic scatter ing of photons or X-rays 
from electrons bound in atoms ( Cooped . Il985l ; IWilliamsl . 
Il977l ). The Compton effect actually played a substan- 
tial role in the early validation of qua ntum ideas, and 
with the experimental work of DuMond (jDuMondl . Il947f ) 
became a practical tool for the investigation of electron 
momentum distributions. Modern experiments involve 
energies in the region of 10 to hundreds of keV, i.e. ener- 
gies that in general are very large as compared to atomic 
Fermi energies (eV). 

The observable in Compton scattering, the so-called 
Compton profile, is the longitudinal momentum distri- 
bution of the initially bound electrons, in direct corre- 
spondence with the scaling function F (y) determined in 
quasi-elastic scattering from nuclei, see Sec. IVIIl The 
momentum distribution of the electrons exhibits a more 
complex structure, as not only the bound-electron n(k) 
play a role, but also — and often more prominently — 
the conduction electrons for metallic targets. 

As an example we show in Fig. [35] the Compton pro- 
file of Sodium. The contribution of the core electrons 
has been removed; the parabola-like part below a mo- 



FIG. 35 Compton profile of Na as a function of electron mo- 
mentum, with the contribution of the core electrons removed 
and finite-resolution effects unfolded (jEisenberger et all . 

dzi). 



When compared to quasi-elastic scattering on nuclei, 
Compton scattering possesses an additional feature: for 
crystalline samples, the Compton profile can be measured 
for many different orientations of the probe. This allows 
one to derive the reciprocal form factor B(r), from which 
one can extract the autocorrelation function in r-space, 
a quantity that helps to deduce the spatial structure of 
the molecules in the crystal. 

A variant of Compton scattering is electron Compton 
scattering, where the photon is replaced by an electron, 
typically in the 50keV energy region. While the produc- 
tion of a good beam and the detection of the scattered 
particle is much easier, the requirement of very thin tar- 
gets, unfortunately, partially offsets these advantages. 

A second area where inclusive scattering represents a 
popular tool is the scattering of low-energy neutrons from 
condensed matter systems such as quantum liquids. Neu- 
trons with energies ranging from lOmeV to lOkeV, from 
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reactors or spallation neutron sources, are employed. The 
energies are again very high as compared to energy scales 
of the system, of order mcV for liquid Helium for exam- 
ple. 

Similar to quasi-elastic scattering from nuclei, the 
measured structure function S(q, lo) yields, in the impulse 
approximation, an integral over the momentum distribu- 
tion. In the limit of large |q|, it can be written as a 
function of one variable Y, the longitudinal momentum 
component. 

In quasi-elastic neutron scattering, emphasis has been 
placed on the understanding of the role of FSI, for two 
reasons. First, the interaction between two atoms is 
fairly singular at short inter-atom distances r; the typi- 
cal Lennard- Jones potential rises very steeply (r -12 ) at 
small r. This strong FSI has a pronounced effect upon the 
observables, to the extent that the longitudinal momen- 
tum distribution can only be extracted after correction 
for FSI. 




FIG. 36 Profile for neutron scattering from liquid Helium at 
0.35K. The curves, calculated using GFMC and including 
FSI effects, are shown for differe nt fractions o f the B ose con- 
densate, 10% giving the best fit l|Sosnik et qZ.I . Il990h . 



FSI was studied in great detail for a second reason. 
Much of the emphasis in the field was placed on the mea- 
surement of the fraction of Bose condensate in supcrfluid 
4 He. This Bose condensate was expected to produce a 
(5-function like spike in the momentum distribution n(k) 
at k = 0, which would lead to a spike in the response at 
Y = 0. This feature was not seen, a fact that now has 
been understood as a consequence of FSI. 

The understanding of FSI in neutron scattering has 
many aspects that are parallel to the discussion given in 
Sec. IIIII In particular, it has been found that the main 
effect of FSI is a folding of the IA response. The folding 
function has a width governed by the atom-atom total 
cross section. It has also been found, that for a quantita- 
tive understanding of FSI, it is imperative to include in 
the description of the initial state the atom-atom corre- 
lation function g(\fl ~ rj|). The treatment of the atom- 



atom potential is comparatively di fficult, yet is p ossible 
with hard-core perturbation theory (jSilverl . Il988l ) . 

As an example, we show in Fig.[3S]the scaling function 
measured at a mome ntum transfer of 23 A on super- 
fluid Helium at 0.35K (|Sosnik et aZ.l . ll990j) . At this tem- 
perature a Bose condensate is predicted to occur. The 
S (Y=0) -function due to the condensate, smeared with the 
FSI folding function, explains the data for a condensate- 
fraction of 10%. 

The last related area we want to address concerns 
deep inelastic scattering (DIS) of GeV electrons or muo ns 
from nuclcons (for a review see e.g. (|Ellis et a/.l . ll996h h 
Here, the energy spectrum of the inclusively scattered 
lepton is used to derive the momentum distribution of 
the quarks bound in the nucleoli. This process provided 
the first direct evidence for the existence of pointlike 
constitue nts of fractional charge, subject to asym ptotic 
freedom (|Gross and Wilczekl . Il973t iPolitzerj . Il973h . and 
still is one of the main sources of information on nucleon 
structure. 
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FIG. 37 Structure function F2 as a function of the momentum 
transfer, for different v alues of the Bjorken scaling variable x 
(|Eidelman et 01.1,120051 ). 



DIS is generally analyzed in terms of the Bjorken scal - 
ing variable x = Q 2 /2mcu ([Bjorken and Paschosl Il969l ) 
appropriate for constituents of negligible rest mass. The 
variable y, which in the limit of small constituent mass 
equals the Nachtmann variable £ (but for a trivial factor) 
is more appropriate for extending scaling to the lower mo- 
mentum transfers ( Benhar et all l2000f ). In terms of £, 
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DIS and quasi-clastic scattering from nuclei have many 
things in common. 

Apart from the parallel aspects (see Sec. |TTJ) DIS ex- 
hibits one special feature: the evolution of the scaling 
function with increasing Q 2 . This evolution is a conse- 
quence of the fact that the lepton, at larger and larger 
Q 2 , resolves more and more of the nucleoli structure 
and the other nucleon constituents, the gluons. This 
evolution of the structure function has been studied in 
great detail, and is well under stood in terms of Quantum 
Chromo Dynamics (QCD) (jAltarelli and Parisil Il977t 
iGribov and Lipatovl . Il972f ). Fig. [37] displays the evolu- 
tion of the proton structure function with momentum 
transfer. 

Unlike e.g. the neutron-scattering mentioned above, 
DIS is always analyzed in IA, neglecting the FSI of the 
recoiling quark; the structure functions then are inter- 
preted directly as the quark distribution functions. The 
effects of FSI, a lthough known to be present even in the 
Q 2 — > oo limit ( Brodskv et all f200^) . are neglected de- 
spite the fact that mo del calculations indicate that the y 
are of substantial size ( Paris and Pandharipandel l2002t ). 

Analogous experiments on DIS of neutrinos have also 
been performed; due to the small rates, the data base in 
this area is much more restricted. 



XIV. CONCLUSIONS 

The field of inclusive quasi-elastic electron-nucleus 
scattering has seen important progress during the last 
decade, both in terms of experimental results and theo- 
retical understanding. 

Experiment has greatly benefited from the high- 
intensity GeV-energy facilities and the high performance 
spectrometers and detectors that became available. 
This has allowed one to extend the data to extreme 
values of momentum transfer and energy loss. As 
a consequence, we now have, at least for selected 
nuclei, a fairly complete data base (see the Web page 
|http : //f acuity . V irginia . edu/ qes- archive that 
gives a rather complete collection of the available cross 
sections). Not yet satisfactory is the situation for the 
longitudinal strength at large momentum transfer, 
where the data base is very narrow. Improvements 
would also be desirable at low momentum transfer, 
where much of the data base comes from experiments 
done in spectrum-acquisition mode. 

Theory has also made considerable progress. Nuclear 
many-body theory today provides reliable spectral func- 
tions, which are at the basis of any quantitative under- 
standing of quasi-elastic scattering, and especially impor- 
tant in the region of high momentum transfer and not 
too large energy loss. In particular, the nucleon-nucleon 
short-range correlations, which have long been known to 
play a major role, are now included in an adequate fash- 
ion. 

The y-scaling analysis of the data clearly shows that 



clastic scattering off a single nucleon is the dominant re- 
action mechanism at x > 1. This is the region where 
quantitative information on nuclear properties can be ex- 
tracted. 

For the treatment of the rccoil-nuclcon final state inter- 
action, various approaches have been developed, applica- 
ble in different regions of momentum transfer and energy 
loss. The results of calculations carried out within the 
scheme widely adopted at large momentum transfer (typ- 
ically q| > 1 GeV), based on the eikonal approximation, 
indicate that final state interaction effects are large, in- 
deed dominant, in the low energy loss tail of the inclusive 
cross section, and their inclusion leads to a quantitative 
account of the existing data up to x ~ 2. 

The role of meson exchange currents, which surpris- 
ingly have a large effect at rather low momentum trans- 
fer and across the entire quasi-elastic peak, also has been 
much better understood. Issues not yet satisfactorily re- 
solved concern the final state interaction at large momen- 
tum transfer and very low energy loss, corresponding to 
x > 2, and the role of non-nuclconic degrees of freedom 
at very large momentum transfer. Recent studies of FSI 
in (e, e'p), car ried out the with i n the formalism described 
in Section [nil (Schi aviIla et aLl . l2005f ). show that the spin 
dependence of the NN scattering amplitude plays an im- 
portant role. The possible relevance of these effects in 
inclusive processes requires further investigations. The 
treatment of knock-out processes accompanied by exci- 
tation of the nucleon also is not yet entirely satisfactory. 

Recently, much consideration has been given to the 
connec tion between electron- and neutrino nucleus scat- 
tering (|Amaro et all 120051 : iBenharj . 120051 ). The gener- 
alization of the existing theoretical approaches and the 
exploitation of the measured (e, e') cross sections to pre- 
dict the analogous neutrino-induced reactions, such as 
(z/, e) or (y, [i), will be of great importance to reduce the 
systematic uncertainty in the interpretation of neutrino 
oscillation experiments. 

Only partially exploited is the relation to quasi-elastic 
scattering of hadronic probes (p, 7r, K, ..) from which ad- 
ditional, and largely complementary, information could 
be learned. 
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